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ABSTRACT 

Context. A survey of mid-infrared gas-phase emission lines of H2, H2O and various atoms toward a sample of 43 embedded low-mass 

young stars in nearby star-forming regions is presented. The sources are selected from the Spitzer "Cores to Disks" (c2d) legacy 

program. 

Aims. The environment of embedded protostars is complex both in its physical structure (envelopes, outflows, jets, protostellar disks) 

and the physical processes (accretion, irradiation by UV and/or X-rays, excitation through slow and fast shocks) which take place. The 

mid-IR spectral range hosts a suite of diagnostic lines which can distinguish them. A key point is to spatially resolve the emission in 

the Spitzer-IRS spectra to separate extended PDR and shock emission from compact source emission associated with the circumstellar 

disk and jets. 

Mefhods. An optimal extraction method is used to separate both spatially unresolved (compact, up to a few hundred AU) and spatially 

resolved (extended, thousand AU or more) emission from the IRS spectra. The results are compared with the c2d disk sample and 

literature PDR and shock models to address the physical nature of the sources. 

Results. Both compact and extended emission features are observed. Warm (rex few hundred K) H2, observed through the pure 

rotational H2 S(0), S(l) and S(2) lines, and [S i] 25 fim emission is observed primarily in the extended component. [S i] is observed 

uniquely toward truly embedded sources and not toward disks. On the other hand hot {T^,, > 700 K) H2, observed primarily through 

the S(4) line, and [Neii] emission is seen mostly in the spatially unresolved component. [Feii] and [Sin] lines are observed in both 

spatial components. Hot H2O emission is found in the spatially unresolved component of some sources. 

Conclusions. The observed emission on >1000 AU scales is characteristic of PDR emission and likely originates in the outflow 

cavities in the remnant envelope created by the stellar wind and jets from the embedded young stars. Weak shocks along the outflow 

wall can also contribute. The compact emission is likely of mixed origin, comprised of optically thick circumstellar disk and/or 

jet/outflow emission from the protostellar object. 

Key words. Stars: formation - Stars: low-mass - Stars: protostars - ISM: jets and outflows - ISM: lines and bands - ISM: photon- 
dominated region (PDR) 

1 . Introduction The envelope surrounding young stars is subjected to highly 

energetic phenomena, in particular X-rays and UV radiation and 
Understanding the nature of deeply embedded protostars is a ma- streams of gas ejected from the star-disk boundary region. Both 
H jor challenge in the study of low-mass star formation. In its ear- the slow and high velocity gas can shape the protostellar en- 
5y , liest stages, a young star is embedded in its natal envelope of yelopes on a variety of spatial scales. The related physical pro- 
cold gas and dust, prohibiting direct detection with optical and cesses, outlined below, have a strong impact on the physical 
near-IR instruments. During the evolution of such protostars the structure and the molecular and atomic contents of the surround- 
envelope is dissipated partly due to the ongoing accretion onto jng material. Mid-infrared spectroscopy, the central focus of this 
the central star and disk and partly due to the action of the out- paper, provides a powerful tool to characterize and quantify the 
flows driven by such stars. This embedded phase lasts only a nature and energetics of these interactions as well as the compo- 
short time, ~ 6 x 10^ yr (JEvans et al.||200g), but is critical for sition and temperature structure of the material involved, 
the subsequent evolution of the star: it is the phase in which the 

final mass of the star, the overall mass of the circumstellar disk Based on previous infrared and submilUmeter data, a general 

(and thus its ability to form planets), and the initial physical and picture has emerged of the structure of e i nbedded young stellar 

chemical structure of the disk are determined. objects (YSOs) (e.g., "Young et al."2004'; 'j0rge nsenet"alll2OO2l 

Chiang et al. 2008; To bin et al. 2010 ). Take as an example 
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IVelusamv et alJ 120071 and references cited) and studied in de- 
tail with submillimeter molecular lines as well as dust contin- 
uum emission (see summary in van Kempen et al. 2009b). The 
young star is surrounded by a roughly spherically symmetric en- 
velope with a mass of ~5 M© and a power-law density structure, 
with n(H2) * 10^ cm"^ at 1000 AU from the star, decreasing to 
~ 10^ cm"-' at -7000 AU. The infrared images reveal prominent 
red- and blue-shifted outflow lobes traversing the envelope with 
an opening angle up to 1 10" and continuing out to 0.2 pc. These 
outflow lobes are produced by jets and/or wide angle winds 
propagating at speeds up to 200 km s"', which interact with the 
surroundings and create a shocked shell of material th at is en- 
traine d and prop elled outward at some speed (e.g.. Ballv et alJ 
I2007t iRav et al.ll2007t IShang et al.|[2007h . The most prominent 
example of this interaction is the bow shock at the tip of the out- 
flow. In addition, blobs of shocked warm gas and dust are seen 
at various locations along the outflow axis, which may be inter- 
nal shocks in the (variable) wind or jet itself. Similar physical 
structures are expected for other embedded YSOs, although the 
densities at 1000 AU are up to an order of magnitude lower for 
sources with lower mass envelopes ( J0rgensen et al. 2002). 

The outflow cavities create a low-extinction pathway 
for ultraviolet (UV) photons produced in the accreting gas 
column at the star-disk boundary layer to escape out to 
much further distances dSpaans et alJ Il995t iBruderer et alJ 
2009b a) than poss ible in a purely spherical symmetry (e.g. 



StauberetaT] l2004 . For a UV luminosity of 0.1 Lq typ- 



ical of low-mass YSOs, the enhancement of the UV ra- 
diation field with respect to the general radiation field is 
about Go ~ lO'* at 100 AU. Moreover, UV radiation is created 
along the outflow by the shock s themselves, especially by 
the higher velocity -/-type shocks dHollenbach & McKedl 19891; 
iNeufeld & Dalgarnolll989h . Evidence for UV enhancements Go 
of a few hundred out to 0. 1 pc have recently been inferred 
from observations of narrow high-7 CO submillimeter e mission 
lines sur rounding the outflow cavity walls (v an Kempen et alJ 
I2009blal see their Fig. 13). Thus, the Spitzer IRS sfit of 
about 5 - 10" will pick up emission from gas with a range 
of densities from 10"* - 10^ cm"-', exposed to UV radiation 
with Go ~ 10^ - 10"*, and subjected to shocks with velocities 
up to several hundred kms"' for typical cloud distances of 
100- 400 pc. 

The consequences of UV radiation impacting on dense gas 
have been studied for decades through models of so-called 
Photon -Dominated Regions (or Ph otoDissociation Regions, 
PDR) ("Tielens & Hollenbachl ll985blli iHollenbach et al.lfl99H 
[Hollenbach & Tielens 1999). In the PDR layers closest to the 
U V source most molecular material is dissociated and the atomic 
species with ionization potential below 13.6 eV are ionized. This 
layer is characterized by bright fine structure lines such as [Oi], 
[Fe i], [Fe ii] and [Si ii]. [S i] is never strong as most of the sulfur 
is in the form of S^ in the warm PDR. At deeper layers, H2 can 
survive, yet remains highly excited, r esulting in strong pure ro- 
tational H2 lines (e.g. Kaufman et al.ll2006, applied to regions of 
massive star- formation). 

Shocks impacting the surrounding envelope can be divided 
into two types: 7 - and C-type shocks. The so-called dis- 
sociative 7-shocks (lpraine&McKedll993HHollenbachlll997t 
IWaknslev et al.ll2005h heat the gas up to 10^ K. Infrared emis- 
sion comes from the dense hot post-shock gas in which most 
of the molecular material is still dissociated and a large frac- 
tion of the atoms are ionized resulting in strong emission 
lines of singly and doubly ionized atoms (e.g., [Fen] 26 and 
35 fim, [Niii] 6.6 fim). This includes atoms which cannot 



be directly photoionized by phot ons with energies <13-6 eV 
such as the [Neii] 12.8 /urn line (IHollenbach & McKed [19891: 
Hollenbach & Gorti 2009). Further downstream, molecular hy- 
drogen forms in the cooling gas and strong high-7 pure rota- 
tional line emission is expected. [Ne 11] emission is strong only 
in high velocity shocks (vs > 70kms"'), while the strength of 
other fine structure lines (e.g. [Si] 25 fim, [Fei] 24 yum and 
[CI i] 1 1.4 ^m) does not strongly d epend on the shock veloc- 
ity (see IHollenbach & McKedl 19891 for details). Note that these 
fast shocks imply travel times through the envelope of only 
(< lO^yrs), unless the material is significantly slowed down by 
the in teraction. For the so-called non-dissociative C- shocks (e.g. 
Neufeld & Kaufman 1993: Kaufman 1995: Kaufman & Neufeldl 



1996h . the temperatures never become high enough to dissociate 
molecules; these shocks show strong lines from excited H2, CO, 
H2O and OH (e.g. Kaufman & Neufeld 1996). 

To date, PDR and shock models have mostly been ap- 
plied to the general interstellar medium (ISM) near massive 
young stars or supernova remnants. Previous observations of 
low mass young stellar objects (YSOs) using ground-based ob- 
servatories or the ISO -SWS instrument (de Graauw et al. 199g 
Ivan den Ancked IT 9991) either lacked the sensitivity to observe 
the emission lines or the spatial resolution to separate the var- 
ious emitting regions (see Ivan Dis hoeck 200 4 for an overview) . 
The sensitive InfraRed Spectrograph (IRS) (Houck et al. '20()4|) 
on board the Spitzer Space Telescope (Wern er et al. 2 004) brings 
the detection of the mid-IR H2 pure rotational and atomic fine 
structure lines within reach for young Sun-like stars in nearby 
star-forming regions. The combination of high sensitivity, mod- 
erate spectral resolution R = A/AA - 600, and modest spatial 
resolution makes Spitzer well suited for the study of the gas in 
the ISM around low-mass young stars in the nearest (< 300 pc) 
clouds. Moreover, the Spitzer aperture allows us, using an op- 
timal extraction method, to disentangle both spatial components 
from the IRS echelle and long-slit spectral images down to scales 
of a few arcsec, better than can be achieved with single-dish mil- 
limeter telescopes. The combination of line fluxes, line ratios 
and spatial extent can help to identify the dominant heating pro- 
cesses - PDR or shock-driven - taking place, and thus the phys- 
ical structure of the envelope on scales of a few hundred to a 
few thousand AU. The Spitzer data thus complement traditional 
studies with (sub)millimeter telescopes. 

Theory predicts that disks are established early in the em- 
bedded phase of star formation, as the inevitable byproduct 
of the collapse of a rotating core (e.g. iTerebev et al.l 11984 
ICassen & Moosmanl[l981i) . Such 'young' disks, some of which 
show Keplerian velocity patterns, ha ve been detected in sub- 
millimeter interferome try data (e.g . Keene & MassonI 119901 ; 
Lommen et al. 2008; J0rgensen et al.|[2009l) . but little is known 
about their physical and chemical structure. Detailed models of 
X-ray and UV irradiation of disk surface layers exist for more 
evolved pre-main seque nce stars where the envelopes have been 
comp l etely dissipat ed (Bergin et al. 2007; Gorti & Hollenbach] 
I2OO8I: IWoiflce et al.l 12009.: .Bergia .2009.) . Spitzer observations 
have revealed surprisingly rich spectra of such protoplanetary 
disks, with the lines seen either in absorption in a near edge- 
on geome try (Lahuis et al. 2006) or in emission (Carr & Najit^ 
l2008; Sa lvk et al.l2008i) . Lines from H2O, HCN, C2H2, CO2 and 
OH have been detected, indicating hot molecular gas at tem- 
peratures of ~ 500 - lOOOK arising from the inner few AU of 
the disk. This hot gas is also seen in ground based CO infrared 
data (e.g. Najita et al. 2003 ; Blake & Boogert 2004 ). Such disks 
show various atomic transitions as well, most prominently the 
[Nen] 12.8;um line dLahuis et al.ll2007bt iPascucci et afl 120061; 
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lEspaillatetani2007h (see also models bv lMeiierink et al.ll2008l) . 
An interesting question is therefore whether young disks in the 
embedded phase studied here have similar spectra as the more 
evolved disks, or whether they show additional features, for ex- 
ampl e due to the accretion shock of material fal ling on to the 
disk jPontoppidan et al.ll2002l; IWatson et al.ll2007l) . 

We present here an overview of the mid-IR gas-phase 
lines detected in embedded sources, including class and 
class I sources observed in the Spitzer legacy program "From 
Molecul ar Cores to Plane t Forming Disks" ("Cores to Disks" 
or c2d) dEvans et al.ll2003h . This program has collected a large 
sample of IRS spectra toward sources in the nearby Chamaeleon, 
Lupus, Perseus, Ophiuchus, and Serpens star-forming regions. 
High-5/A^ spectra have been obtained within the 5 - 38 //m 
range for 226 sources at all phases of star and planet formation 
up to ages of ~5Myr, of which 61 are presented here. Many 
of the sources are also included in the s eries of papers dis- 
cussing the ice features in the c2d sample dBoogert et al.ll20()8t 
iPontoppidan et al.ll2008t fOberg et al.ll2008h . In Sections Hand [3] 
the source selection and data reduction are explained. In Section 
|4]the observed atomic fine-structure and H2 emission lines and 
the derived parameters are presented. Section previews the re- 
sults within the context of PDR and shock heating and by com- 
paring the embedded sources in the present sample with the more 
evolved disk sources in the larger c2d program that have been 
discussed in tLahuis et al..i2007bi) . 

2. Observations 

The data presented in this study were selected from the sample 
of IRS spectra observed within the Spitzer c2d legacy program. 
The c2d IRS component consists of two programs of compara- 
ble size, referred to as the first- and second-look programs. The 
first-look program (PID #172) was restricted primarily to known 
low-mass young stars, embedded YSOs and pre-main-sequence 
stars with disks with stellar masses M <2 Mq and ages < 5 Myr, 
and a sample of background stars. A few Herbig Ae stars were 
included as well. The c2d source selection criteria were defined 
to be complementary to those of the Spitzer legacy program 
"The Formation and Evolution of P lanetary Systems" (FEPS, 
iMever & the entire FEPS teaml2007h . The second-look program 
(PID #179) was, for the most part, devoted to IRS follow-up 
spectroscopy of sources discovered in the IRAC and MIPS map- 
ping survey s, including a newly discovered cluster of young stars 
in Serpens dHarvev et al.ll2006l) . For all first-look observations, 
the integration times for the short-high (SH) and long-high (LH) 
modules (R = 600, 10-37/im) were chosen such that theoretical 
signal to noise ratios (S /Ns) of at least 100 and 50 on the contin- 
uum were obtained for sources brighter and fainter than 500 mJy, 
respectively. Deeper integrations were not feasible within the 
c2d program. Spectra taken using the short-low (SL) or long-low 
(LL) modules {R = 60 - 120, 5-14 //m and 14-38 fim respec- 
tively) always reach theoretical S /N ratios greater than 100. For 
the second-look IRS targets similar S /N limits were obtained 
wherever possible. However, since the second-look contained a 
number of very weak sources (down to a few mJy) this was not 
always achieved. 

The 61 sources presented in this paper were all selected 
from those showing the IQfim silicate band in absorption. 
Thi s criterion includes em bedded Stage I sources (as defined 
by iRobitaille et al.l l2006l) plus obscured disks (e.g. edge-on 
stage II sources (i > 65 degrees), Crapsi et al. 2008), such 
as CRBR 2422.8-3 423 dPontoppidan et al.ri200l and IRS 46 
dLahuis et al.ll2006 '). All of the sources have a rising SED in the 



mid-IR, the characteristic for Class I YSOs. The selected sources 
are listed in Table[T]which gives the basic observation and source 
parameters plus the adopted distances. Note that in most cases, 
it is not possible to determine whether the infrared source is 
dominated by an envelope or a disk without additional spa- 
tially resolved far-infrared and/or (s ub)millimeter data . A sub- 
sequent submillimeter line study by van Kem pen et al.l d2009ch 
of Ophiuchus L1688 has been able to separate truly embedded 
YSOs (Stage I) from obscured disks or confused sources. Of the 
61 sources, 43 are found to be truly embedded; 10 sources are 
obscured stage II sources (edge-on disks or disks behind large 
columns of gas and dust); and for 8 sources the classification is 
unclear. Some obscured disks wer e not included in the c2d IRS 
survey of gas-phase lines in disks ( dLahuis et alj2007bl) . so their 
data are presented here for completeness, but are analyzed sepa- 
rately. 

3. Data reduction 



The c2d reduction pipeline (Chapter 3 of lLahuisll2007l and the 
c2d-IRS delivery documentatiorQ) was used to reduce the IRS 
data, starting from the S15, S16 and S17 archive data. The same 
c2d pipeline products as included in the final c2d Legacy data 
deliver}0 were used for the spectral line analysis (see Sec. 13.4b . 
Most of the analysis focused on the SH and LH data, since for 
intrinsically narrow lines the SL and LL data are generally lim- 
ited by the low line/continuum ratio. The SL data were included 
in the analysis to search for higher-excitation H2 v = lines, in 
particular the S(3) transition. 

If a module was not included in the c2d observation the 
Spitzer archive was searched for additional pointed observations 
to supplement the missing data. For some sources data are miss- 
ing, in particular SL. For example, Serp-SMM3 and Serp-SSM4 
have not been observed in the low resolution staring mode (see 
TablelTJ and therefore do not cover the higher rotational {J > 3) 
H2 lines. Some of these are covered through mapping observa- 
tions, however these are not included in this a nalysis. 

The c2d disk sources dLahuis et al.l2007bh are used for com- 
parison with the embedded sources and have all been repro- 
cessed using the same redu ction to provide the d ata presented in 
this paper. Compared with Lahuis et al.l (|2007b) this has led to 
an increased number of detection for the disk source s, most im- 
portantly new detections of [Ne n] dGudeletal.l2010h and [Fe n]. 

3.1. Separating disk and cloud emission - optimal extraction 

A major concern when studying YSOs in the mid-IR is the possi- 
ble contribution of extended (envelope or local cloud), yet struc- 
tured, emission in the sometimes complex star-forming regions. 
Variations in the spatial distribution of the emission (see Figure 
[TJ, both in the continuum and in spectral lines, limit the use of 
'sky' observations to correct for extended emission components. 
For this reason the c2d team developed an optimal extraction 
algorithm for IRS pointed observations starting from the SSC 
Archive droop products. 

The Spitzer sUt size is ~ 4 - 5" for the SH module (10- 
19.5;um) and ~ 7 - 10" for the LH module (19-37;um). At 
a distance of 100 pc this corresponds to physical sizes of 
~ 400 - 500 AU and ~ 700 - 1000 AU, respectively. The clouds 



' http://data.spitzer.caltech.edu/popular/c2d/ 
2006 1 20 1 _enhanced_v l/Documents/C2D_Spec tr_ExpLSupp.pdf 

- The c2d legacy data are accessible at 
http://ssc.spitzer.caltech.edu/legacy/c2dhistory.html. 
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Fig. 1. Spitzer-IRS SH minimap observed around VSSGl in Ophiuchus. The left plot displays a Spitzer-IRAC image at 4.5 //m (including the 
high excitation H2 S(9), S(10) and S(ll) lines) showing the star plus diffuse extended emission. The right plot shows the spectrum observed in 
the on-source Spitzer-IRS observation (top panel) and the off-source spectra obtained with the Spitzer-IRS minimap. In gray the total observed 
emission and in black the unresolved source emission toward VSSGl are presented. Virtually all H2 S(l) 17.03 pm and S(2) 12.28 pm emission is 
extended in this source. [Neii] 12.81 pm emission is absent toward VSSGl. 



observed in the c2d program are assumed to be located at dis- 
tances ranging from 125 pc (Ophiuchus) to 260 pc (Serpens) 
increasing the physical area observed. The full IRS aperture in 
the SH and LH spatial direction is ~ 2.5 - 3 times larger than 
the beam size. At the observed cloud distances this means that 
the aperture probes physical scales of several thousand AU. This 
makes it ideally suited for detecting warm or shocked H2 emis- 
sion from the extended (remnant) envelope, outflows, or the dif- 
fuse local cloud emission. For the sources studied in this work 
however the local cloud emission will potentially confuse com- 
pact disk emission. Distinguishing between compact (disk and 
inner envelope) and extended (remnant envelope, outflow, or dif- 
fuse cloud) emission is therefore of vital importance for study- 
ing the emission lines originating in the circumstellar disks or 
the cores of jets/winds launched on small scales. 

The optimal extraction developed by the c2d team uses an 
analytical cross-dispersion point spread function (PSF) to de- 
scribe the source profile (see Figure |2]i. The PSF plus an ex- 
tended emission component are fit to the observed crosstalk or 
straylight corrected echelle (SH and LH) and longslit (SL and 
LL) images. The PSF is described by a sine function with a har- 
monic distortion component (resulting in a significant broaden- 
ing of the profile wings as shown in Figure |2]i. The wavelength 
dependence of the PSF parameters (the order trace, the width, 
and the harmonic distortion) are characterized using a suite of 
high-5/A^ calibrator stars. The extended emission component is 
approximated by a low order polynomial folded with the flatfield 
cross-dispersion profile. The flux calibration is derived from 
the calibrator stars using Cohen templates and MARCS models 
(iDecin et al.ll2004l) provided through the Spitzer Science Center. 
More details about the characterization and cali bration of the c2d 
optimal extraction is presented in Chapter 3 of Lahuis (2001) and 
in [Lahuis et al. (2007b). The optimal extraction returns the total 
flux (the source flux plus extended emission in the IRS beam) 
along with an estimate of the extended emission component with 



error estimates. The S /N of the extended emission component 
can vary significantly depending on the quality of the raw im- 
age data and on deviations of the extended emission from the 
assumed uniformity across the IRS slit. 

Source sizes as listed in Table |2] are determined from 
the width of the observed cross dispersion source profile (see 
Figure |2]i deconvolved with the average width of the standard 
star source profiles. 



3.2. SH mini maps 

In an early phase of the c2d project H2 and [Ne 11] lines were 
detected toward some of the c2d sources. Initial results using 
an early version of the c2d optimal extraction, from IRS start- 
ing observations, suggested some of the observed emission to be 
extended. Therefore, as part of the second-look program, follow- 
up mini-maps were taken using the SH module to check for ex- 
tended emission at offsets positions of ~10-15" with respect to 
the sources. Five maps were defined to include off-source obser- 
vations around eig ht sources. F igure [T] shows the results toward 
VSSGl, whereas .Lahuis et al.1 (12007b.) (their Figure 1) show the 
results toward the disk sources Sz 102 (Krautter's star), Haro 1- 
17, and EC 74. Here, Figure[T]presents the observed H2 0-0 S(l), 
S(2), and [Ne 11] emission from the first-look on-source observa- 
tions together with the off-source observations from the second- 
look mini maps. It is seen that most of the H2 emission is ex- 
tended and that it may vary on scales of a few ", especially 
for the S(l) line. However, some emission lines such as [Ne ii] 
are typ ically seen to be limited to the source itself (Lahuis et al.l 
I2007bh . 

The mini map results confirmed the results from the c2d op- 
timal extraction procedure. It showed that it was possible to re- 
liable extract the extended emission component from the IRS 
staring observations. The procedure was further developed and 
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Fig. 2. Illustration of the Spitzer IRS cross dispersion profile used in the optimal extraction (see Sec. l3.lt . The left plot shows a fit to the data of 
an IRS SH order of a pointsource with a moderate but clear sky component in the IRS spectra. The SH pixel size is ~2.3 arcsec. For LH the profile 
is similar but the pixel size is ~2 times larger. The observed data are plotted with black plusses. Overplotted are the total profile fit (source plus 
extended emission) and the included extended component. The shape of the extended emission reflects the IRS flatfield of the (for this source) 
assumed uniform extended emission. The right plot shows a comparison of an IRS PSF profile (Sine plus harmonics) for a pointsource compared 
to the profiles of an undistorted Sine, a Gaussian, and a Lorentzian profile with the same FWHM. Note the significant variation in the strength and 
shape of the profile wings. The correct characterization of both the width and the wings of the profile for all IRS orders is essential for extracting 
the proper source and sky spectra. 



has since been used in several studi es, e.g . lLahuis et alj (l2007bh : 
iGeers et alJ (l2006l) : lBoersma et alJ (12008 



3.3. 1-D spectra 

After extraction, the 1-D spectra are corre cted for instrumental 
fringe residuals (" Lahuis & Booger3l2003h . order matching (by 
a few %) is applied, and a pointing flux-loss correction is per- 
formed onthe compact source component (see Chapter 3 of 
lLahuisll2007h . Pointing offsets up to a few arcsec can have a no- 
ticeable impact on the derived fluxes of lines observed with the 
SH and SL modules, e.g. H2 0-0S(1), S(2), [Neii], and [Neiii]. 
For example, dispersion offsets within the nominal 3cr pointing 
uncertainty of Spitzer (~1" for medium accuracy peakup) can 
lead to SL and SH wavelength dependent flux losses up to ~ 10%. 
Fortunately for all targets a combination of modules is available 
which permits us to apply a reliable pointing flux loss correction. 



unresolved component. Spectral lines are fit to obtain line fluxes 
for both spatial components. Uncertainty estimates, as listed in 
Tables [3] and |4] are derived from the residuals after line fitting, 
or, in the absence of a spectral line, using the line width derived 
from the instrumental resolution. An uncertainty estimate is de- 
rived for both the extended and compact line fit and combined. 
As a result, the 1-cr uncertainty estimates can vary widely for 
sources with a similar continuum flux and integration time. This 
may, for example, be the result of the presence of artifacts result- 
ing from hot pixels or small variations in the extended emission 
which are not fully accounted for in the spectral extraction. 

Typical mean 4cr uncertainties prior to subtraction of 
the extended component for the high resolution modules 
range from ~0. 1-2x10"''' ergcm^^s"' with extremes of 
~5xlO"'^ erg cm"- s"' and ~lxlO"''^ erg cm"^ s"' . Note that the 
uncertainties used in this paper include an extended component 
uncertainty. 



3.4. Spectral analysis 

The IRS SH and LH modules cover the positions of the three 
lowest pure rotational lines of H2, fine structure emission from 
[Neil] (12.8 //m), [Neiii] (15.55 //m), [Fei] (24//m), [Feii] (17.9 
and26.0jum), [Si] (25.25 jum), [Siii] (18.7 //m), [Sin] (34.9, fim), 
and hundreds of H2O transitions at a resolving power of R - 
A/AA = 600. The IRS SL module covers the higher rotation 
H2 lines up to 7/ = 7 at a resolving power of 60 - 120. Line 
fitting and flux integration is done to all H2 and atomic lines 
using routines from osi^j. 

The extended spectral component (either the estimate or a 
Gaussian fit depending on the final S /N in the corrected spec- 
trum) is subtracted from the total spectrum to isolate the spatially 



^ osiA is a joint development of the ISO-SWS consortium. 
Contributing institutes are SRON, MPE, KUL and the ESA 
Astrophysics Division, http://sws.ster.kuleuven.ac.be/osia/ 



4. Results 

A search for the major atomic fine-structure lines and H2 lines 
in the Spitzer-IRS wavelength range has been performed for all 
61 sources in our sample. Line detections are shown in Figures 
|3]and|4|for the embedded sources and described in Sections |4T| 
and l4.2l In Figure|5]the measured line strengths and upper limits, 
converted to solar luminosities, are shown as function of Lik (the 
continuum luminosity between A - 12.8 - 15.5 i-im), see Section 
14.41 Figure |6| shows the ratios and ratio limits for a number of 
lines. Lines of molecules H2O, C2H2, HCN and CO2 have also 
been searched for and are discussed in Section l43] 

For both the atomic lines and those of H2 the observed upper 
limits are in most cases of the same magnitude as the median line 
fluxes of the detections (see Figure |5]l. It is therefore likely that 
for most sources emission in one or more lines will be present 
but has gone undetected. The upward trend observed in the upper 
limits above an IR luminosity of ~ lO^^L© is the result of the 
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Fig. 3. Detections of [Ne II], [Fe II], [S I] and [Si II] at the > Act level toward the c2d sample of embedded YSOs. Overplotted in red is a Gaussian 
fit to the observed spectrum. Sources indicated with a star show evidence of a (partial) contribution of extended emission to the total line flux. 



dominance of target shot noise in the IRS detector noise at higher lections ~30% show emission from more than one line. Emission 
source signals. of most atomic lines is observed as both compact and extended 

emission. 

4.1. Atomic fine-structure iines 

4.2. Moiecular iiydrogen 
Figure|3]shows all 4cr detections of atomic lines. One or more of 

the [Neil], [Feu] 18 and 26yum, [Si], and [Sin] emission lines Figure |4] shows all 4cr H2 emission line detections toward the 

are observed toward ~2/3 of the 43 embedded sources. The mea- 43 embedded sources in our sample. Emission from one or more 

sured line fluxes of the extended and compact components and H2 lines observable with the Spitzer-IRS (S(0) to S(7)), has been 

the upper limits are listed in Table[3] Of the sources with line de- detected toward 36 out of 43 sources, with 3 1 out of 43 showing 



Fred Lahuis et al.: c2d IRS spectra of embedded sources: gas-phase emission lines 



Hj S(0) Hj S(l) Hj S(8) H3 S(3) H, S(4) H, S(5) H, S(6) H^ S(7) 



Hj S(0) Hj S(l) Hj S(2) Hj S(3) H3 S(4) H, S(5) H3 8(6) H, S(7) 



450 -tSS I 



.200 , ..150 



fv'^/vv 



J r 



1500 ^20 I at 



'25 I ■ . 



.40 .25 I X 



.400 4.20 I * 1; 



350 I ..300 .'150 

yM h i * "yv~^ vyv^ioi 



.40 I . .7 I X, 



.35 . .45 



•lfft\^ ' 



.aoo ..15 



■iTM^ 



r^ 



.20 I M . 



*«M/1 




.70 , ■..ao , ■(..,0 , m 




.20 , . .5 



'■^ 



,200 h K .■300 ',40 I 



n,A 



.350 ..50 . 45 



150 I .'60 .,70 

iJlig \.yv\fli **JM^ 

' 1 l M | 1 1 — I— I 1 — I — I— 



..7 I ..35 



.8000 ..1000 ..700 



■vi*sv)" ijwv "iry!* 



■fc««-u>\ 



25 I X 



.40 I «.'" I « 

■ - t . . ■ ill A 







-f^^/V 



4-JL_. 



■yVfwM 






ilA/y^V 



v4^ 



ftfjiyfiy 



Wv/'lw 









ijAiyA 



■nHUjV 



wyvv 



ir^vM 



H£: 



vV'vV 



Jft^ 



Hf Try nr tT" 



ifei 






Ik 






M^ 



3 



1*. 



Jffi 



o^v^ 



ite 



iife 



iT-- 



m 



350 H 

diJLii 



lyy*- 



iJVV 



riiXli^ 



n^ 



-aka 



lArAw 



|J.„„ Juui 



<Yw* 



■WV 



ftiMf 



irt>%#i 



«iA|/J< 



ify^Af 



vv'iWV 



(iV^/^v 



Liui 



<ii|fy%* 



jftAj* 



lyVNJV 



ilr<fMJb 



iW^Y 



<fin\pi\f- 



Yv^ 



-H 1 1— 



■f/V^ 






H 1 1— 



i^ 



J*s 



*t^ 



°L 



jm 



Jv 



Mt 



r-M 



H 1 h- 






:X 



350 1 X..30 I X.''^ I X' 

,JL- J — JL 



350 I 

mi 




.400 I * . 



VKfA^ 



350 I X . 

juJLlj. 



\4|A>#>* 



0/**— ' 



/Vy 



vW^M 




150 ,. 



-■^rfyr^ 



\>y^ 



40 I * ■ 



(WAV 





so 1 



20 I 



45 I * . .45 I * . .30 1 K . 

AjL# — -JL— . . ■ JL. < 

I " I — I — I — I— I — I — I — I— 



r-W 



kfc 



350 



vr-v/t-y* 



jVV^ 



<V'^*I' 



■>nlT 



.-vJ^W 




««— 'Vw 



f-<./Si.Ai/ 



45 I * . 



"-illjw 



Ip-f-^ 



f^f^ 



rJw*l 



VD^iMl 



■yVT- 



(l-JljS-i 



*JV/f«V 



■yyv^ 



4IAMi 



fY7 



rMf^ 



ma^ 



1 )K 

hMdld 
TP-Inn 



vMU 



-Ji-vV 



150 



S/q|l~«.. 



150 |L 



VMU 



Ifs/l" 






, J|l]Uf 



l/iVM 



tuiMrtl 



ill^ 



ji 



ft 



Ac 



Jfc. 






200 



yv-^ 



■yVy^ 



■Ai«Art 



■^.A>~« 



UWu 



Wfrt** 



WN/lr 



lA-VV 



—I — I — I— 



t^/\^ 



jyw; 



w\/* 



iWi* 



A-VSI 



irr 



ifsM/V 



Jfc 



27.7 26.7 16.8 17.2 12.1 12.5 8.9 10.5 7.2 



7.4 5.6 6.6 5.0 6.0 



27.7 26.7 16.8 17.2 12.1 12.5 8.9 10.5 7.2 8.8 6.4 7.4 5.6 6.6 5.0 6.0 



Wavelength [^m] 



Wavelength [/im] 



Fig. 4. Detections of H2 emission lines at the > 4cr level toward the c2d sample of embedded YSOs. Overplotted in red is a Gaussian fit to the 
observed spectrum. Sources indicated with a star show evidence of a (partial) contribution of extended emission to the total line flux. 
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low-y lines and 19 out of 43 high-/ lines. The H2 line fluxes 
for the compact and extended components and the line flux upper 
limits are listed in Table |4] The warm gas traced by the low- J 
lines is observed primarily in the extended component while the 
hot gas probed by the high-/ lines is observed mostly in the 
compact component. 

The warm gas (rex< 700 K) traced by the lower rotational 
lines (Ji < 2) is mostly extended, except for a few cases. Hot 
iTex~ 1000 - 1500 K) gas is traced by the higher rotational lines 
(Ji > 4). For 7/ > 5 no estimate of extended emission is available 
because the IRS is spatially undersampled at these wavelengths 
and the usability of the optimal extraction is limited. Mapping 
observations (e.g. Neufeld et al. 2009) show that, in general, the 
H2 high-y lines have a similar distribution to the lower excitation 
lines. Since most of the 7/ = 4 detections show no strong signs of 
an extended emission component, we assume all hot H2 emission 
to have a compact origin. 

In the simplest analysis, the H2 exc itation is assum ed to be 
in local thermal equilibrium (LTE) fe.g. jThi et al.ll200l1) with an 
ortho-to-para ratio determined by the kinetic temperature of the 
gas (following Sternberg & Neufeld 1999). For gas temperatures 
100, 150, and >200K, the ortho-to-para ratios are 1.6, 2.5, and 
3, respectively. Assuming optically thin emission, the integrated 
flux of a rotational line 7,, —> Ji for a given temperature Tex is 



Fui(T,.) = p-N(n2)AuiXu(T,,)n erg s"' cm-^, (1) 

where A is the wavelength of the transition, A^(H2) the total col- 
umn density, A„; the spontaneous transition probability, and Q 
the source size. For high enough densities (n > lO^cm"^), the 
population Xu follows the Boltzmann law 



Xii\i ex/ — 



gN{2Ju + l)e 



-Ej/kT,, 



Q(T,,) 



(2) 



where Ej is the energy of the upper level, ^n is the nuclear sta- 
tistical weight (1 for para and 3 for ortho H2), and QiTg^) the 
partition function for the given excitation temperature Tex. 

Using the above equations, excitation temperatures, column 
densities and H2 gas masses, or limits on these, can be de- 
rived from the observed line fluxes and upper limits. A source 
size has to be known or assumed to derive column densities. 
Observations of multiple H2 lines in PDR s and shocks typically 
show multiple temperat ure components (Ivan den Ancker et alJ 
ll999HHabart et al.ll2005h . e.g. a warm (rex of a few 100 K) and 
a hot (rex ~ lOOOK or higher) component. In our data, often 
only a few lines per source are detected. If 7/ < 2 lines are de- 
tected, a (limit of the) temperature of the warm gas is derived. 
For a few sources more than one of these lines are detected in 
the extended component. The derived temperatures vary from 
85 to 120 K with uncertainties of ~ 20 K. Therefore we adopt a 
temperature of lOOK for the warm component for all sources. If 
two or more higher excitation lines (7/ > 2 and higher) are de- 
tected a temperature for the hot component is derived; otherwise 
a temperature of lOOOK is adopted. This is done for both spatial 
components. For some sources multiple H2 lines are detected al- 
lowing a temperature estimate to be made. Temperatures for the 
hot component, varying from ~ 500 - 1500 K, are observed (see 
Figure |7]i in both spatially extended and compact emission. In 
the analysis no extinction correction is applied and only H2 is 
included in the mass calculation. 

The column density averaged over the IRS aperture is de- 
rived using above equations. For all compact source emission 



the source size is smaller than the IRS aperture (Sec. [3]). Since 
this is unknown, a typical size r = 50 AU is adopted (gas masses 
up to a few lO"-' Jupiter masses, see below and Figure |8] give 
an upper limit to the emitting area of 50 AU assuming a den- 
sity n > 10^ cm""*). The fitted or adopted excitation tempera- 
ture plus the 7up column densities provide the total column den- 
sity for that temperature, which in turn gives the total H2 gas 
mass, Mh2 -nr- xN x 2mii/Mj with ma - 1 .674 ■ 10"^"^ gr and 
Mj = 1.9 • 10^° gr, excluding a mass contribution by He. Note 
that the derived column density scales as 1 /r^ while the derived 
gas mass is independent of the adopted beam and/or source size. 
Therefore the derived column densities can only be used as order 
of magnitude estimates and are listed for both the warm and hot 
gas in Table |5] Figure [8] shows the H2 mass distribution of the 
warm (Tex - 100 K) and hot (Tex = 1000 K) components for the 
truly embedded sources. H2 masses of a few Jovian masses for 
the warm gas and a few milli- Jovian masses for the hot gas are 
observed. Hatched bars include all sources where H2 is detected. 

4.3. Molecular emission 

Features of H2O, C2H2, HCN and CO2 emission are searched for 
in the full sample. Detecting these species is difficult because 
the emission features are generally weak and ride on top of a 
strong mid-IR continuum. H2O emission is most easily detected 
as the spectra contain a large number of both spectrally unre- 
solved and blended emission lines across the complete IRS-SH 
spectral range. Thus, a full spectral range correlation can be ap- 
plied. Detecting C2H2, HCN and CO2 is more difficult because 
only the Q-branches can be observed (P- and R-branch lines are 
undetectable at the IRS resolution). The Q-branches are, more- 
over, weak and partially blended with H2O lines, as seen for 
example in the spectra of VSSG17 and L1689IRS5 spectra in 
Figure|9]where the C2H2, HCN and CO2 Q-branch positions are 
marked. 

Toward 12/43 truly embedded sources spatially unresolved 
H2O emission is positively identified and the presence of C2H2, 
HCN and CO2 is, by visual inspection, observed toward a small 
number of sources (see, e.g., the spectrum of VSSG17 in Figure 
|9]l. Toward 4 sources H2O is clearly detected in both the opti- 
mal extraction and the full aperture extraction spectra. For 11 
more sources H2O is positively detected in the optimal extrac- 
tion spectrum and tentatively detected in the full aperture spec- 
trum. Figure |9] shows the continuum subtracted SH spectra with 
a LTE water model overplotted for two of the embedded sources. 
The observed emission features are similar to those seen in th e 
more evolved disks around TTauri stars bv ISalvk et al.l (l2008l) : 
ICarr&Naiital (12008ft . 

To derive basic excitation parameters the full IRS SH spec- 
tra are compared with single te mperature UTE models us- 
ing the HITRAN 2004 database (' Rofliman et all l2005h . The 
model assumes a single slab of gas of solid angle (Q), tem- 
perature (T), column density (A^) and intrinsic linewidth (b). 
Without exception the derived excitation temperatures are high, 
Tex = 700 ~ 1500 K. Best fits are achieved for models with large 
H2O column densities, A^ > 10'^ cm"^, and small emission re- 
gions, r < few AU, comparab l e to w h at is found for the TT auri 
stars observed by Salvk et al.1 (l2008h : ICarr & Naiital (l2008l) . An 
accurate and fully constrained fit is not possible with the low 
{R = 600) resolution SH spectra. It may be po ssible to extract 
more information through non-LTE models (see lMeiierink et al] 
2009) but such models are outside the scope of this paper. 

For C2H2, HCN and CO2 none of the spectra have sufficient 
S IN to derive reliable excitation temperatures and column den- 
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Fig. 5. Line luminosities and upper limits (in solar luminosity) as 
functions of the mid-IR luminosity (calculated from the IRS spec- 
tra) for the estimated extended and compact source contributions 
to the total line strength. The red circles show the tru ly embedded 
source s, white stars the disk sources (this paper and lLahuis et al.l 
l2007bh and blue squares the unclassified sources. Shown are all 
4cr line detections, the symbol size scales with line error. The ob- 
served upper limits listed in Tables|3]and|4]are shown as grayscale 
images with dashed contours. The upper limits are presented as 
number density distributions with the contours bracketing the 10 
and 90 % ranges. All sources plotted left of the vertical dashed 
line either have no detected continuum in the IRS spectra or the 
spectrum does not cover the mid-IR range (x-position offset for 
clarity). 



Fig. 6. Line ratios (filled symbols) and line ratio upper and lower 
limits of [Ne ii] with a number of H2 and atomic lines. The use of 
symbols and colors is the same as in Figure|5] All sources plotted 
left of the vertical dashed line either have no detected continuum 
in the IRS spectra or the spectrum does not cover the mid-IR range 
(x-position offset for clarity). The two horizontal dotted lines span 
the range of 0.001 to 10 within which most detections fall. This 
range is also indicated in Figure[TO]which shows [Ne 11] line ratios 
from shock models by Hollenbach & IVIcKee ( 1989). Note that 
most line ratios diplayed here are upper and lower limits. Only 
for [Ne ii]/H2 S(4) a reasonable number of line ratios (filled red 
circles) exist. 



sities or even solid detections. High spatial and/or spectral res- 
olution observations are required to further constrain the model 
p arameters. The im portance of high resolution data is illustrated 
in lSalvk et al.l ( 12008,) who obtained high resolution 3/im spectra 
for two disk sources, AS 205 A and DR Tau, following the initial 
detection of hot H2O emission in the Spitzer-IRS SH spectra. 
We therefore do not list fit parameters for individual sources but 
limit ourselves to summarizing detections in Table |2] 



4.4. Correlations and line ratios 

Figure |5] shows the observed line luminosities (for all 4cr detec- 
tions) as function of the mid-IR (12.8 - 15.5 //m) continuum lu- 
minosities as calculated from the Spitzer-IRS spectra. The mid- 
IR continuum luminosity is used as a qualitative measure of the 
hot-dust component. 

Both the line strengths and the upper limits fall in a narrow 
band, ~2 orders of magnitude, spanning the full Lir range. All 
species show a positive trend with Lir, which is most clearly 
observed in the compact H2S(4) and [Neii] line emission. As 
discussed above, the presence of weaker line emission than cur- 
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Fig. 7. Excitation diagrams of the observed H2 emission. Plotted 
with black circles is the observed spatially compact emission 
and with red squares the spatially extended emission. The ar- 
rows are the derived upper limits. Overplotted with respectively 
dashed and dashed/dotted lines is a two temperature fit to the 
observed H2 emission, as seen by the difference in slope going 
from the low to the high energy states. The inferred temperatures 
are listed in Table |5] Note that the low temperature component 
is not always well constrained. 



rently observed cannot be excluded. However, the current dataset 
does exclude strong line emission for the weak IR sources. 

Figure |6] shows the line ratios and line ratio limits of [Neii] 
with a number of H2 and atomic lines. All ratios fall within a 
range of a few orders of magnitudes and show no obvious trends. 
The line ratio limits do not reveal any sources deviating from 
this, consistent with the fact that the observed line luminosities 
are close to the line luminosity upper limits as shown in Figure 
|5] Higher resolution and high sensitivity observations (e.g. with 




Fig. 8. Distribution of H2 masses (in Jupiter mass) for both the compact 
and the extended component assuming T^^ = 100 and 1000 K. The 
gray bars represent upper limits on the gas mass, while the hatched 
bars indicate sources where H2 is detected at 4cr or more. For 100 K 
this includes sources with the H? S(0) and/or S(l) line detected and for 
1000 K sources with any of the higher J lines detected. The highest bar 
also includes all sources with masses higher than the upper plot limit. 



0.07 - 




12 14 16 18 

Wavelength [/7,m] 

Fig. 9. Molecular emission toward two truly embedded sources. Plotted 
for each source are the continuum subtracted Spitzer IRS-SH spectra 
(top spectrum in black) and a single temperature LTE water model at 
Tev = 900 K (bottom spectrum in red). Vertical lines indicate the posi- 
tions of a few important spectral features; from left to right H2 S(2), 
[Neil], C2H2, HCN, CO2 and H2S(1). In VSSG17 CO2 and tentative 
C2H2 and HCN emission can be seen. 



the MIRI instrument aboard the James Webb Space Telescope) 
will be required to overcome this limitation. 
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5. Discussion 

5.1. General trends 

Atomic fine structure lines, H2 pure rotational lines and molec- 
ular emission features are observed toward the majority of the 
sources in our sample, both as compact and extended emission. 
Table |2] presents a qualitative overview of the identified emis- 
sion. The compact emission is dominated by [Neii], high-/ H2 
emission lines and hot H2O emission. The extended emission 
is dominated by low-7 H2 emission lines. A small number of 
sources show line emission from both the extended and com- 
pact components at the same time. H2O emission is seen in the 
compact component toward both compact and extended emis- 
sion dominated sources. 

The increase of the line strength with increasing mid-IR lu- 
minosity suggests that the heating mechanism of the dust and the 
excitation of the line may have the same physical origin. Besides 
heating due to the protostellar luminosity, UV radiation from the 
star-disk accretion boundary layer escaping through the outflow 
cones may heat both the dust and the gas (see §5.3). 

A global comparison in Figures [5] and [6] between the trul 



embedded sources and the disk sample of (ILahuis et 



•n tne truly 
^ I2007hh 



shows that both samples exhibit spatially compact and extended 
emission at the same time. Both samples also show compact hot 
water emission. The main difference is in the statistics of the ex- 
tended emission. While in the compact emission the difference 
between the two samples is small (e.g. hot water emission is ob- 
served toward ~40 % of the disk sources and ~30 % of the truly 
embedded sources) the difference is large for the extended com- 
ponent. Extended emission is observed toward ~l/3 of the disk 
sources and ~2/3 of the truly embedded sources. The extended 
emission observed toward the disk sources represents the PDR 
emission from the local environment, in particular the extended 
PDR in p Oph. Therefore it is fair to assume that a significant 
fraction of the extended emission observed toward the truly em- 
bedded sources has its origin from the source and not from the 
environment. 

Because the Spitzer-IRS spectral and spatial resolution is 
moderate, the level to which we can qualitatively and quan- 
titatively constrain the excitation and the physical environ- 
ment of young stars is limited. High spatial and/or high spec- 
tral resolution observations are required to separate the dif- 
ferent emission components either through direct observations 
or by spectral decomposition. A prime example of this is the 
study of [Ne 11] emission toward stage II disk sour ces. The ini- 
tial Spitz er detections of compact [Ne 11] emission dLahuis et alJ 
[2007b; Pascucci et al.l !2007) raised excitement and suggested 
a possible disk surface origin of the [Nenl emission whic h 
was supported by model results from iGlassgold et al.l (l2007h . 
However, recent high resolution ground based observations 
dHercze g et"ani2007l;lPascucci & Sterzi5l2009l:lvan Boekel et alj 
[2009) have shown that high resolution foll ow-up observations 
are es sentia l for a proper i nterpretation. iPascucci & SterziH 
(120091) and Ivan Boekel et alj (12009.) argue for an outflow/jet 
origin of the spatially unresolved Spitzer [Ne 11] emission ob- 
served toward TTau, Szl02 and VWCha. Recent model results 
(Hollenbach & Gorti 2009) confirm that observable [Ne 11] emis- 
sion may be produced through X-ray and FUV irradiation of the 
proto-planetary disk surface and through i nternal shocks in pro- 
tostellar winds. This is confirmed by Gii del et alj (|2010) where 
disk and jet contributions are separated through a statistical anal- 
ysis of [Ne It] emission of a large sample of Spitzer sources. 

Although a detailed analysis of individual sources may 
not be possible, general conclusions may be drawn. By 



looking at the overall line fluxes and line ratios observed 
(see Figure s|5] and |6]l we can draw general conclusions about the 
protostellar environment. In the following sections we will there- 
fore compare the observations with published shock and PDR 
models and discuss possible origins of the observed emission. 

5.2. Shock excitation 

Many of the observed emission lines, including [Neii], [Fen], 
[S i] and [Si 11] are predicted to be strong in 7-type shocks 
(Hollenbach & McKee 1989) with the main exception being 
the low-y H2 lines. The hot (~1500K) H2 emission observed, 
in particular, is expected i n the cooling post-shock gas (e.g. 
Ivan den Ancker et alJ Il999h . The observed strengths of the 
atomic emission lines are comparable with predictions by shock 
models, but the observed line ratios are not sufficiently con- 
strained to draw any conclusions about precise shock veloci- 
ties. However for [Neii] to be strong, shock velocities in excess 
of ~ 70kms"' are required. Figure [TOl shows that the range of 
line ratios predicted by the models is consistent with those ob- 
served (Figure ^ over a range of expected pre-shock densities 
and shock velocities. 

The H2 S(0) line is predict ed to b e the weakest /-shock fea- 
ture in Hollenbach & McKe^ (Il989t) . The observed line ratios 
and upper limits of the [Ne 11I/H 2 S(0) line rat io for the extended 
component in the disk sources dLahuis et al.ll2007h) is <1 (see 
Figure |6]l. In /-shock models the [Nen]/H2S(0) ratio is up to 
two orders of magnitude higher (see Figure [TOl l. We therefore 
exclude /-shocks as the origin for the extended low-7 warm H2 
emission toward the embedded and disk sources. The extended 
[Fell] emission observed toward the disk sources may have a 
similar origin as the warm H2 emission though this cannot be 
constrained by the /-shock models itself. 

In non-dissociativ e C-shocks low-/ H2 emission is strong 
(iKaufman & Neufeldl fl996) and the observed strengths of the 
H2 lines are consistent with C-shock model predictions. At the 
same time atomic lines are not expected to be strong. The low-/ 
H2 lines, mostly observed in the spatially extended component, 
may partly have a C-shock origin in material entrained by the 
outflow. However given that [Fe 11] is also observed as extended 
emission, C-shocks are not likely to be the sole origin of the 
extended line emission. 



5.3. PDR models 

Models of photodissociation regions (PDRs) predict line 
strengths for the low-/ H2, [Fen] and [Sin] emission compa- 
rable to the observed line strengths fo r n ~ 10^ - lO^cm"^ 
and a rad iatiori field of ~ 10' - 10^ Go dHoUenbachet al.ll 19911: 
Kaufman et al.ll2006l) . These values are consistent with those es- 
timated for the outer envelope. The H2 excitation temperatures 
of ~ 350 - 700 K (see Figure|7|l are comparable with gas surface 
temperatures for PDRs with the range of dens ities and radiation 
fields mentioned above (lKaufmanetal.ll 19991 and the PDRlQ). 

A PDR origin of the observed extended emission therefore 
appears to be plausible. It is likely that the observed emission 
originates in the PDR tracing the surfaces of the extended out- 
flow cavity for the embedded sources, although there may also 
be a contribution from extended lower density cloud material 
in the beam for some sources, especially in Ophiuchus. Weak 
C-shocks may be present as well (see Section 15.2b . The pres- 



"^ Photo Dissociation Region Toolbox 
http : //dustem . astro . umd . edu/pdrt/ | 
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Fig. 10. [Nell] line ratios from y-shock models bv lHoUenbach & McKed ([T98p) for pre-shock densities of 10' (black) and IC* (red). Note that 
for the low density case the [Ne ii] ratios with the low-7 H2 lines are large and quickly fall off the plot. The hatched area indicates the range of 
observed line ratios (see Figure|6ll. 



ence of extended photon-heated gas along outflow cavities has 
recently also been inferred from narrow high-7 CO submil- 
limeter emission lines ( van Kempen et al. 2009bllatlSpaans et al.l 
Il995h and has also bee n suggested ba^sed o n ISO-LWS atomic 
fine-structure hne data (iMolinari et al.ll2001h . Although the [S i] 
emission is almost exclusively observed in the extended com- 
ponent, [S i] is not strong in PDR models. At the warm PDR 
surfaces most sulfur will be in the form of S^. 

The origin of the exciting UV radiation field is still un- 
der debate (see §1). Suggested origins are UV photons: i) es- 
caping fromjflie_star-disk boundary layer near the central pro- 
tostar (ISpaans et al.l 119951) . and ii) from high-velocity shocks 
directly associated with the more extended je ts and outflows 
dNeufeld & Dalgarnoll98"9tlMohnari et al. '2001). In both cases, 
the UV field consists of a continuous spectrum with possibly 
some discrete lines due to recombining atoms superposed. In our 
sample we observe sources simultaneously exhibiting signatures 
of extended PDR/C-shock and compact /-shock emission and 
sources showing only extended PDR/C-shock emission and no 
compact shock emission (see Table |2]i, thus leaving all options 
open. 

In order to advance on the current situation it will be im- 
portant to get multiple, when possible spatially and/or spectrally 
resolved, detections per source on key lines such as [Ne ii], [S i], 
[Oi] 63.2yum, [Cii] 157.7 /zm and the low-/ H2 lines. [Neii] and 
[S i] are clear indicators of 7-shocks, whereas the [O i]/[C 11] ra- 
tio and the absence or presence of [Si 11] and [Fe 11] help to distin- 
guish and_characterize PDR and C-shock emission. In his PhD 
thesis lvan den Anckeil(ll999L Chapter 10) presents a useful sum- 
mary of the relevant issues. High velocity resolution [Neii] ob- 
servations are possible from the ground although the lines pre- 
sented here are not strong enough for current instrumentation. 
JWST will provide an inventory of the key mid-IR emission lines 
towards a large sample of embedded sources. Herschel will pro- 
vide the possibility to observe highly-excited CO, [O i] and [C 11] 
and provide high velocity resolution observations for some far- 
IR lines through HIFI. 



5.4. Molecular emission and absorption 

While molecular emission is seen toward a number of sources, 
none of the truly embedded sources show absorption bands of 
HCN, C2H2 and/or CO2. This in contrast with hi gh and interme- 
diate mass stars JLa huis & van Dishoeck 2000; iBoonman et al.l 
120031: iBoonman & van Dishoeckl |2003l Thi et al. in prep.) 
where these molecules are characteristic of h ot-core chemistry 
(iDotv et al.ll20"0l iRodgers & Charnlevll2003h . Toward ULIRGs 
abs orption of C2H2 , HCN and CO2 has been observed by 
Lah uis et aP (l2007ah suggesting the presence of pressure con- 
fined massive star formation. The absence of signatures of C2H2 
and HCN absorption towards the low-mass sources in our sam- 
ple is not surprising, however To produce absorption features 
observable with IRS one needs both high column densities, 
requiring orders of magnitude enhancement of the molecular 
abundances compared with the normal cloud abundances, and a 
strong infrared background. These conditions are met in the in- 
ner envelope around higher-m ass stars during the hot-core phase 
(iLahuis & van Dishoeckl2000l) . but for solar-mass stars the inner 
envelope regions do not become hot enough (T > 300 K) to pro- 
duce the required chemical enrichment. Only in the inner regions 
of protoplanetary disk are these conditions met and absorp- 
tion features have only be en observed with Spitzer and near-IR 
ground-based instruments (; Lahuis et aT] l2006': 'Gibb et al.ll2007h 
under a favorable, near edge-on, source alignment toward two 
sources (IRS 46 and GV Tau). 

Molecular emission from H2O is observed in a quarter of the 
truly embedded sources in our sample and emission from C2H2, 
HCN and CO2 is tentatively observed in a few sources. What is 
the origin of this emission: the young disk or outflow? 

Mid-infrared H2O and OH emission is observed to be asso- 
ciated with outflows such as HH 211, where prominent lines are 
detected at the bow sh ock position offset from the protostar it- 
self (Tap pe et al.ll2008l e.g.). Emission features of C2H2, HCN 
and CO2 seen toward Cepheus A East are evidence of gas-phase 
abundance enhancemen ts through sputterin g of icy grain man- 
tles in shocked regions (ISonnentrucker et a l. 2006, 2007). Since 
our data show neither the high-7 OH lines nor prominent C2H2, 
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HCN or CO2 emission, these scenarios do not appear to explain 
our data. 

Strong H2O mid-infrared emission has also been observed 
with Spitzer towar d the stage source (NGC1333-IRAS4B by 
IWatson et al.ll200'^ and attributed to disk accretion shocks, but it 
is unclear whether this interp retation is unique and how genera l 
this phenomenon is (see also lj0rgensen & van Dishoeckll2010l) . 
None of the other stage sources in our sample (in partic- 
ular VLA.1623-243, IRAS_16293-2422, IRAS_16293-2422B, 
Serp-S68N, Serp-SMM3 and Serp-SMM4) show hot H2O emis- 
sion. One Stage I source, GSS30-IRS1, exhibits stro ng CO mid- 
infra red lines indicative of accretion shocks ( Pontoppid an et al] 
12002 ) but does not show any H2O lines in our Spitzer data. 

On the other hand, sti'ong H2O, C2H2, HCN, CO2 and/or 
OH lines have clearly been detected from a number of more 
evolved stage II sources arising from the inner AU of the disk 
(ICarr & Naiita'2008':'Salvk et al.'2008). The molecular emission 
observed toward the embedded sources in our sample has char- 
acteristics similar to those observed toward these disk sources: 
high Tex, high n and a small size of the emitting region. In par- 
ticular, Tex is significantly higher tha n the value of ~ 170 K ob- 
served toward NGC 1333 IRAS4B dWatson et alJl2007h . It is 
therefore plausible that the molecular emission observed in our 
sample also has its origin in the inner AU of the young disks. 
In other words, we may be looking through a window onto the 
embedded young disks (see e.g. ,Cernicharo et aL i2000i). 



Envelope 




Fig. 11. Cartoon image to illustrate the various emission components 
for the embedded sources as seen by mid-IR Spitzer-\KS observations. 
Extended emission may be observed as PDR and/or C-shock emission 
from the irradiated and accelerated gas in the cavity and envelope sur- 
faces. Compact /-shock emission can be observed from shocks within 
high- velocity jets, indicated by filled circles. Disk emission may be ob- 
served for systems viewed along the outflow cavity. 



5.5. Source environment 

The general picture which emerges is that of a young star sur- 
rounded by a forming disk embedded in the natal envelope. 
Winds and jets from the young star produce holes and cavities 
in the envelope. Within the jets and in its boundary regions with 
the envelope, supersonic shock waves are produced which lo- 
cally compress and heat the gas producing, among others, hot 



H2 and [Ne 11] emission. Non-dissociative C-shocks in material 
entrained by the extended outflow produce warm spatially ex- 
tended H2 emission. C-shocks are however not the sole origin 
of the observed extended emission: the presence of spatially ex- 
tended [Fe 11] and [Si 11] emission argues against this and points 
to the presence of extended PDRs. UV radiation produced by 
the young star system heat the cavity walls and envelope surface 
creating extended PDR regions. Emission from the young disk, 
in particular that from molecules such as H2O, can be observed 
directly by looking along the outflow cavities and/or through a 
patchy envelope. 

Spatially extended PDR and C-shock emission is observed 
toward ~2/3 of the sources, including the stage sources. Some 
of this emission may come from the local environment, as ob- 
served toward the disk sources. Extended emission is observed 
toward ~l/3 of the disk sources. Therefore a significant fraction 
of the extended emission observed toward the embedded sources 
may originate from the sources themselves. This is consistent 
with most sources having a sizeable outflow cavity. The outflow 
cavities are most easily observed toward sources where the out- 
flow is seen face-on or at moderate inclination angles. Toward 
systems seen at large inclination angles the extinction through 
the envelope will prohibit detecting the emission from the out- 
flow cavity walls, though this will depend on the observed wave- 
length of the emission and its location within the cavity emis- 
sion. Thus, given a uniform line-of-sight distribution, extended 
emission from the outflow cavity walls should be observed to- 
wards a large fraction of the sources. This is the case in our sam- 
ple of embedded sources. 

The presence of hot H2O emission suggests a strong simi- 
larity between the embedded stage-I disks and the more evolved 
stage-II disks. None of the known stage-0 sources in our sample 
show signs of hot H2O emission, implying that their disks have 
either different characteristics or are much more obscured. 

Figure[TT]shows a graphical representation of the picture de- 
scribed above displaying the different components and the loca- 
tion of the observed emission. In this picture, the main difference 
with the more evolved stage II disks is that the latter sources lack 
the extended emission from PDR and C-shocks along the cavity 
edge. The limitations of the Spitzer-IRS data (low spatial and 
spectral resolution) prohibit us from making conclusive state- 
ments on individual sources. 



6. Conclusions 

A survey of the mid-infrared gas phase pure rotational lines 
of molecular hydrogen and a number of atomic fine structure 
transitions has been carried out toward a sample of 61 class I 
sources with the Spitzer-IRS of which 43 are truly embedded 
protostars. Both extended and compact source emission has been 
detected toward a majority of the sources. Compact emission, 
most prominently seen in hot H2, [Ne 11] and H2O lines, is seen 
toward ~l/4 to 1/3 of the truly embedded sources. The compact 
hot H2 and [Ne 11] emission is believed to have its origin predom- 
inantly in dissociative 7-shocks in the outflow jets, though a true 
disk component may be present. The H2O emission is consistent 
with an origin in the inner regions of the young protostellar disk, 
similar to the H2O emission from the more evolved systems. 
Extended emission, dominated by low-7 H2 and [Feii] emis- 
sion, is observed toward ~2/3 of the truly embedded sources, 
both stage and I sources. The extended emission is thought to 
originate from PDRs along the outflow cavity walls with possi- 
ble contributions from non-dissociative C-shocks. 
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These observations have shown the capability of Spitzer to 
observe the gas around young stars. It clearly shows emission 
from PDR and shock excited gas by the presence of various H2 
and atomic lines. However, its limited spatial and spectral resolu- 
tion prevents a quantitative analysis of the excitation conditions. 

Sources with atomic fine structure line emission, H2 and 
H2O emission are excellent candidates for follow up with var- 
ious ground based instruments and near-IR and far-IR space- 
based instruments on Herschel and JWST. The major far-IR 
cooUng Unes of [Oi] and [Cii] as well as high-/ CO observ- 
able with Herschel can provide important excitation constraints 
(Ivan Kempen et al." 2010). Several of the sources in this sample 
are targets in the "Water in Star-forming regions with Herschel" 
(WISH) and "Dust, Ice and Gas in Time" (DIGIT) Herschel key 
programs. Future mid-IR observations with JWST-MIRI will al- 
low the detection of weak emission lines and provide many more 
diagnostics to separate the physical processes. Such high spec- 
tral and spatial resolution observations are essential to constrain 
the location and extent of the emission (either by direct detection 
or through velocity resolved spectra) to properly understand the 
physical environment around the young stars. 

The spectra of the compact component of the truly embedded 
sources show remarkable similarities with those of the evolved 
disk sources, both in terms of absolute fluxes and flux ratios. 
Both characteristic 7-shock emission, hot H2 and [Neii], and 
true disk emission through the molecular emission lines are ob- 
served. In particular, the molecular emission may offer a unique 
tool to study the young forming disks in their earliest phases and 
compare these with the more evolved stage II disks. 
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Table 1. Source list 



Embedded sources 



AOR 



(SL/SH/LH)*ndith. 



[pc] 



1 


LDN 1448 IRS 1 


3''25"'09\4 


30° 46' 21' 


.7 


5656832 


(14*1/31*2/60*1)*2 


250 


2 


LDN1448NA 


3'- 25- 36^5 


30° 45' 21' 


.2 


5828096 


(60*1/121*2/60*1)*2 


250 


3 


IRAS 03245+3002 


jh 27m 39^0 


30° 12' 59' 


.4 


6368000 


(14*4/31*4/60*1)*2 


250 


4 


L1455 SMMl 


jh 27m 43s 2 


30° 12' 28' 


.8 


15917056 


(60*2/121*2/-)*2 


250 


5 


L1455 IRS3 


3" 28" 00^4 


30° 08' 01' 


.3 


15917568 


(14*4/121*2/-)*2 


250 


6 


IRAS 03254+3050 


3''28"'34\5 


31° 00' 51' 


.1 


11827200 


(14*2/31*2/60*1)*2 


250 


7 


IRAS 03271+3013 


3''30"M5\2 


30° 23' 48' 


.8 


5634304 


(14*2/31*4/60*1)*2 


250 


8 


IRAS 03301+3111 


31,33,112^8 


31° 21' 24' 


.1 


5634560 


(14*l/6*2/14*2)*2 


250 


9 


Bl-a 


31,330, igs 7 


31° 07' 55' 


2 


15918080 


(14*4/31*2/-)*2 


250 


10 


Bl-c 


31, 33m 17s 9 


31° 09' 31' 


.0 


13460480 


(14*2/-/-)*2 


250 


11 


SSTc2dI033327.3+3I0710 


3h33,„27.<.3 


31° 07' 10' 


.2 


15918336 


(60*3/121*2/-)*2 


250 


12 


HH211-mm 


3'' 43'" 56*. 8 


32° 00' 50' 


.4 


5826304 


(-/31*2/60*2)*2 


250 


13 


IRAS 03439+3233 


31,470,05^4 


32° 43' 08' 


.4 


5635072 


(14*1/121*2/60*1)*2 


250 


14 


IRAS 03445+3242 


31,47,141^5 


32° 51' 43' 


.9 


5635328 


(14*l/6*2/14*l)*2 


250 


15 


IRAS 08242-5050 


gi,25m43sg 


-51° 00' 35' 


.6 


5638912 


(14*l/6*2/14*l)*2 


460 


16 


IRAS 08261-5100 


8''27"'38'.9 


-51° 10' 37' 


.2 


5638912 


(14*l/6*2/14*l)*2 


400 


17 


Ced 1 10 IRS4 


11'' 06"" 46^6 


-77° 22' 32' 


.5 


5639680/12692224 


(-/121*2/60*2)*2 


178 


18 


Ced 1 10 IRS6 


lli,07"'09'.2 


-77° 23' 04' 


2 


5639680/12696832 


(-/121*2/60*2)*2 


178 


19 


IRAS 12553-7651 


12''59"'06'.6 


-77° 07' 40' 


.1 


9830912 


(14*1/31*1/60*1)*2 


178 


20 


IRAS 15398-3359 


151, 43m 02^3 


-34° 09' 06' 


.8 


5828864 


(14*2/31*6/60*1)*2 


100 


21 


GSS30-IRS1 


16''26■"21^4 


-24° 23' 04' 


2 


5647616/12699392 


(-/31*1/14*2)*2 


125 


22 


GY23 


16''26"24M 


-24° 24' 48' 


.2 


5647616/12699136 


(-/31*1/14*2)*2 


125 


23 


VLA 1623-243 


16'■26"■26^4 


-24° 24' 30' 


.2 


9828096/12695040 


(-/121*4/60*4)*2 


125 


24 


WL12 


16'" 26"" 44^2 


-24° 34' 48' 


.4 


5647616/12699136 


(-/31*1/14*2)*2 


125 


25 


IRS37 


lgh27mi7Sg 


-24° 28' 56' 


.6 


5647616/12698880 


(-/31*1/14*2)*2 


125 


26 


WL6 


16''27"■21^8 


-24° 29' 53' 


.2 


5647616/12698880 


(-/31*1/14*2)*2 


125 


27 


Elias32 


16''27■"28^4 


-24° 27' 21' 


2 


12664320/12698624 


(-/31*2/60*2)*2 


125 


28 


VSSG17 


Igh 27m 30^2 


-24° 27' 43' 


.6 


5647616/12698624 


(-/31*1/14*2)*2 


125 


29 


IRS63 


16''31"'35'.7 


-24° 01' 29' 


.6 


9827840/12676608 


(-/31*1/14*1)*2 


125 


30 


L1689-IRS5 


16''31"52M 


-24° 56' 15' 


.4 


12664064/12704256 


(-/31*1/60*1)*2 


125 


31 


IRAS 16293-2422B 


16''32"22\6 


-24° 28' 32' 


.2 


15735808 


(60*l/121*4/60*4)*2 


125 


32 


IRAS 16293-2422 


16''32"■22^9 


-24° 28' 36' 


.1 


11826944 


(60*l/121*4/60*4)*2 


125 


33 


RN0 91 


16''34'»29».3 


-15° 47' 01' 


.3 


5650432 


(14*1/31*1/14*2)*2 


140 


34 


SSTc2dJ182916.2+01822 


18''29'"16».2 


0°18'22' 


.7 


13210112 


(14*l/31*2/-)*2 


260 


35 


Serp-S68N 


18''29"'48M 


1°16'42' 


.6 


9828608 


(-/121*4/60*4)*2 


260 


36 


Seip-SMM4 


18" 29"' 56'. 6 


1°13'15' 


2 


9828608 


(-/121*4/60*4)*2 


260 


37 


EC88 


18" 29"" 57'. 6 


1°13'00' 


.5 


9407232 


(14*l/121*l/60*2)*2 


260 


38 


Serp-SMM3 


18" 29"" 59'. 2 


1°14'00' 


2 


9828608 


(-/121*4/60*4)*2 


260 


39 


R CrA IRS5 


19"01"48'.0 


-36° 57' 21' 


.6 


9835264 


(14*l/6*2/-)*2 


130 


40 


CrA IRS7 A 


19" or 55'. 3 


-36° 57' 22' 


.0 


9835008 


(14*l/31*3/60*2)*2 


170 


41 


Ci-A IRS7 B 


19"01"'56'.4 


-36° 57' 28' 


.1 


9835008 


(14*l/31*3/60*2)*2 


170 


42 


Ci-A IRAS32 


19"02"'58'.7 


-37° 07' 34' 


.7 


9832192 


(60*l/31*8/-)*2 


170 


43 


IRAS 23238+7401 


23" 25"" 46'. 7 


74° 17' 37' 


.3 


9833728 


(14*l/31*8/60*2)*2 


250 



Unknown sources 



# Source 



RA 



Dec 



AOR 



(SL/SH/LH)*ndith. 



[pc] 



44 


IRAS 13546-3941 


13"57"'38'.9 


-39° 56' 00' 


2 


5642752 


(14*1/31*2/60*1)*2 


630 


45 


WL19 


16"27""11'.7 


-24° 38' 32' 


.3 


9829888/12692992 


(-/31*2/14*2)*2 


125 


46 


IRAS 16544-1604 


16"57""19'.7 


-16° 09' 23' 


.8 


5826816 


(-/-/14*4)*2 


140 


47 


SSTc2dJ182813. 2+00313 


18"28"'13'.2 


0° 03' 13' 


.0 


13210368 


(14*l/31*2/-)*2 


260 


48 


SSTc2d J 182849.4+00604 


18" 28"" 49'. 4 


-0° 06' 04' 


.7 


13210624 


(14*l/31*l/-)*2 


260 


49 


SSTc2dJ182901. 8+02954 


18" 29"' 01'. 8 


0° 29' 54' 


.2 


13210112 


(14*l/31*2/-)*2 


260 


50 


SSTc2dJ182914.8+00424 


18" 29"" 14'. 8 


-0° 04' 23' 


.9 


13210112 


(14*l/31*2/-)*2 


260 


51 


EC69 


18"29"54'.4 


1°15'01' 


.8 


9407232 


(14*l/121*l/60*2)*2 


260 



Disk sources 



AOR 



(SL/SH/LH)*ndith. 



[pc] 



52 


VSSGl 


53 


IRS 14 


54 


OphE-MM3 


55 


GY224 


56 


WL20S 


57 


IRS42 


58 


CRBR 2422.8-3423 


59 


IRS46 


60 


IRS51 


61 


EC82 



16" 26'" 18'. 9 
16"26'"31'.0 
16" 27" 05'. 9 
16"27"'11'.2 
16"27"'15'.6 
16"27"'21'.5 
16"27"'24'.6 
16"27"'29'.4 
16" 27" 39'. 8 
18" 29" 56'. 9 



-24° 28' 19". 6 
-24° 31' 05". 2 
-24° 37' 08". 
-24° 40' 46". 6 
-24°38'45".6 
-24° 41' 43". 1 
-24° 41' 03". 1 
-24° 39' 16".2 
-24°43'15".2 
1''14'46".7 



5647616/12699392 

12664576 

6370816 
9829888/12698880 

9829888 
5647616/12699648 

9346048 
9829888/12698624 
9829888/12699648 

9407232 



(-/31*1/14*2)*2 

(60*2/3 l*l/-)*2 

(60*5/31*8/60*2)*2 

(-/31*2/14*2)*2 

(-/31*2/14*2)*2 

(-/31*1/14*2)*2 

(14*1/121*1/241*2)*2 

(-/31*2/14*2)*2 

(-/31*2/14*2)*2 

(14*l/121* l/60*2)*2 

i III (200 pc'Comeron 200851 



125 
125 
125 
125 
125 
125 
125 
125 
125 
260 



"Assumed cloud distances; Chamaeleon (178 pc Whittet et al. 1997), Lupus I, IV (150 pc) and Lupus 1 
Ophiuchus (125nc^d^eu^^LJ989[), Perseus (250 pc, see discussion in lEnoch et aU2006l) . Taurus-Auriga (160pc lKenvon et alJ1994l) , 
Serpens (260Dc lStraizvs et alJ1996l) 
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Table 2. Source charateristics 



Embedded sources 



# 


Source 


Warm H2 


HotHj 


[Neil] 


[Fen] 


[Sin] 


HzO" 


outflow'' 


source size 






ext. 


src 


src 


ext. 


ext. 






["] 


1 


LDN1448 IRSl 








y 






y 


- 


2 


LDN1448 NA 


y 


y 


y 


y 






y 


- 


3 


IRAS 03245+3002 


y 


y 










y 


3.1(1.5) 


4 


L1455 SMMl 


y 


y 






- 




n 


0.2(0.2) 


5 


L1455 1RS3 


y 








- 


y 


n 


- 


6 


IRAS 03254+3050 














y 


1.1(0.8) 


7 


IRAS 03271+3013 


y 


y 




y 






y 


0.2(0.2) 


8 


IRAS 03301+3111 






y 








n 


0.3(0.2) 


9 


Bl-a 


y 


y 


y 




- 


y 


y 


- 


10 


Bl-c 




y 






- 




y 


- 


11 


SSTc2d 1033327.3 


y 


y 


y 




- 




y 


0.3(0.2) 


12 


HH211-mm 


y 


- 




y 


y 




y 


ext 


13 


IRAS 03439+3233 














n 


- 


14 


IRAS 03445+3242 








y 


y 




y 


0.2(0.2) 


15 


IRAS 08242-5050 






y 


y 








0.2(0.2) 


16 


IRAS 08261-5100 


y 












- 


0.5(0.4) 


17 


Ced 1 10 IRS4 


y 


y 


y 


y 






y 


- 


18 


CedlI0IRS6 












y 




0.6(0.5) 


19 


IRAS 12553-7651 














n 


- 


20 


IRAS 15398-3359 




y 






y 




y 


0.5(0.3) 


21 


GSS30-IRS1 




y 










n 


0.4(0.3) 


22 


GY23 


y 










y 


n 


0.6(0.5) 


23 


VLA 1623-243 


y 








y 




y 


7.5(2.8) 


24 


WL12 




y 










n 


0.1(0.1) 


25 


IRS37 


y 




y 








n 


- 


26 


WL6 


y 


y 








y 


n 


0.1(0.1) 


27 


Elias32 


y 












n 


- 


28 


VSSG17 


y 










yy 


n 


- 


29 


IRS63 












y 


n 


0.2(0.2) 


30 


L1689-IRS5 












yy 


n 


0.2(0.2) 


31 


IRAS 16293-2422B 


y 






y 






n 


ext 


32 


IRAS 16293-2422 


y 












n 


ext 


33 


RN0 91 






y 






y 


- 


0.4(0.3) 


34 


SSTc2d 1182916.2 










- 


y 


y 


- 


35 


Serp-S68N 


y 


- 






y 




c 


- 


36 


Serp-SMM4 


y 


- 




y 


y 




c 


ext 


37 


EC88 












y 


c 


1.2(0.9) 


38 


Serp-SMM3 


y 


- 




y 


y 




c 


ext 


39 


R CrA IRS5 






y 




- 


y 


c 


- 


40 


CrA IRS7 A 




y 


y 








c 


0.3(0.2) 


41 


CrA IRS7 B 


y 


y 


y 








c 


0.8(0.6) 


42 


CrA IRAS32 




y 


y 


y 


- 




c 


0.5(0.3) 


43 


IRAS 23238+7401 






y 








c 


0.9(0.7) 








Unknown sources 










# 


Source 


Warm H2 


HotHz 


[Neil] 


[Fen] 


[Sin] 


H20» 


outflow"" 


source size 






ext. 


src 


src 


ext. 


ext. 






["] 


44 


IRAS 13546-3941 














- 


0.6(0.3) 


45 


WL19 


y 










y 


n 


- 


46 


IRAS 16544-1604 




- 


- 








- 


- 


47 


SSTc2d 1182813.2 






y 




- 




n 


0.2(0.2) 


48 


SSTc2d 1182849.4 










- 




n 


- 


49 


SSTc2dJ182901.8 










- 


yy 


n 


- 


50 


SSTc2d 1182914.8 










- 




n 


- 


51 


EC69 


y 






y 


y 




c 


ext 








Disk sources 










# 


Source 


Warm H2 


HotH, 


[Neil] 


[Fen] 


[Sin] 


HjO" 


outflow"" 


source size 






ext. 


src 


src 


ext. 


ext. 






["] 


52 


VSSGl 


y 


y 






y 


yy 


n 


- 


53 


IRS 14 


y 


y 






- 




n 


7.7(0.8) 


54 


OphE-MM3 


y 












n 


7.1(1.5) 


55 


GY224 


y 


y 








yy 


n 


- 


56 


WL20S 




- 


y 


y 






n 


0.6(0.2) 


57 


IRS42 














n 


- 


58 


CRBR 2422.8-3423 


y 


y 










n 


0.1(0.1) 


59 


IRS46 














n 


0.2(0.2) 


60 


IRS51 














- 


0.3(0.2) 


61 


EC82 














c 


0.4(0.3) 



''One yindicates a detection in one extraction; two yyindicates sources with detection in both the optimal extraction and the full aperture extraction. 
''Evidence for the presence of outflows based on a visual inspection of Spitz.er-lKhC images; yes, no or confused. 
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Table 3. Observed linefluxes and Icr uncertainties (10 '^ergcm ^s ') 



Embedded sources 



# 


Source 




[Ne 


11] 




[Ne 


Ill] 


[Fe 


I] 24.04 




src'' 


ext" 


upp" 


src 


ext 


upp 


src ext 


upp 


1 


LDN 1448 IRS 1 


_ 


_ 


(30/23) 


_ 


_ 


(43/7) 


_ 


(110/44) 


2 


LDNI448NA 


240 


23 


(20/3) 


- 


- 


(40/28) 


- 


(140/99) 


3 


IRAS 03245+3002 


- 


- 


(6/5) 


- 


- 


(20/20) 


- 


(110/60) 


4 


L1455 SMMl 


- 


- 


(10/8) 


- 


- 


(9/2) 






5 


LI455 IRS3 


- 


- 


(18/2) 


- 


- 


(15/2) 






6 


IRAS 03254+3050 


- 


- 


(24/13) 


- 


- 


(31/4) 


- 


(94/83) 


7 


IRAS 03271+3013 


- 


- 


(19/11) 


- 


- 


(43/29) 


- 


(72/53) 


8 


IRAS 03301+3111 


390 


- 


(76/60) 


- 


- 


(87/76) 


- 


(170/69) 


9 


Bl-a 


150 


- 


(23/4) 


- 


- 


(30/26) 






10 


Bl-c 


- 


- 


(19/17) 


- 


- 


(20/16) 






11 


SSTc2d J033327.3 


870 


- 


(22/11) 


52 


- 


(11/1) 






12 


HH211-mm 


- 


- 


(6/6) 


- 


- 


(4/4) 


- 


(11/10) 


13 


IRAS 03439+3233 


- 


- 


(24/3) 


- 


- 


(9/1) 


- 


(110/53) 


14 


IRAS 03445+3242 


- 


- 


(110/89) 


- 


- 


(90/72) 


- 


(150/87) 


15 


IRAS 08242-5050 


- 


- 


(110/83) 


- 


- 


(75/9) 


- 


(200/90) 


16 


IRAS 08261-5100 


- 


- 


(130/96) 


- 


- 


(76/46) 


- 


(120/50) 


17 


Ced 1 10 IRS4 


180 


- 


(23/2) 


- 


- 


(6/1) 


- 


(35/24) 


18 


Ced 1 10 IRS6 


- 


- 


(22/2) 


- 


- 


(21/15) 


- 


(73/35) 


19 


IRAS 12553-7651 


- 


- 


(31/24) 


- 


- 


(56/8) 


- 


(130/60) 


20 


IRAS 15398-3359 


- 


- 


(22/18) 


- 


- 


(46/10) 


- 


(160/76) 


21 


GSS30-IRS1 


- 


- 


(640/330) 


- 


- 


(1100/160) 


- 


(2400/1300) 


22 


GY23 


- 


- 


(110/68) 


- 


- 


(91/22) 


- 


(140/80) 


23 


VLA 1623-243 


- 


- 


(5/1) 


- 


- 


(2/2) 


- 


(52/48) 


24 


WL12 


- 


- 


(130/53) 


- 


- 


(200/170) 


- 


(330/92) 


25 


IRS37 


290 


- 


(31/25) 


- 


- 


(38/35) 


- 


(110/84) 


26 


WL6 


- 


- 


(78/42) 


- 


- 


(110/88) 


- 


(210/100) 


27 


Elias32 


- 


- 


(43/5) 


- 


- 


(38/32) 


- 


(270/19) 


28 


VSSG17 


- 


- 


(63/31) 


- 


- 


(79/12) 


- 


(170/130) 


29 


IRS63 


- 


- 


(36/28) 


- 


- 


(59/8) 


- 


(230/120) 


30 


L1689-IRS5 


- 


- 


(89/57) 


- 


- 


(120/16) 


- 


(200/100) 


31 


IRAS 16293-2422B 


- 


- 


(4/4) 


- 


- 


(0/0) 


- 


(35/35) 


32 


IRAS 16293-2422 


- 


- 


(3/3) 


- 


- 


(1/1) 


150 


(32/33) 


33 


RN0 91 


710 


280 


(83/59) 


- 


- 


(93/9) 


- 


(210/110) 


34 


SSTc2dJ182916.2 


- 


- 


(37/34) 


- 


- 


(31/5) 






35 


Seip-S68N 


- 


- 


(14/2) 


- 


- 


(6/1) 


- 


(26/23) 


36 


Seip-SMM4 


- 


- 


(3/3) 


- 


- 


(1/1) 


- 


(5/5) 


37 


EC88 


- 


- 


(22/3) 


- 


- 


(42/7) 


- 


(90/58) 


38 


Serp-SMM3 


- 


11 


(4/1) 


- 


- 


(3/3) 


- 


(23/23) 


39 


R CiA IRS5 


1700 


110 


(140/27) 


- 


- 


(160/22) 






40 


CiA IRS7 A 


3700 


- 


(120/69) 


- 


- 


(420/39) 


- 


(1500/680) 


41 


CrA IRS7 B 


960 


71 


(45/15) 


- 


- 


(100/10) 


- 


(560/380) 


42 


CrA IRAS32 


200 


- 


(18/3) 


- 


- 


(23/20) 






43 


IRAS 23238+7401 


150 


- 


(22/12) 


- 


- 


(54/12) 


- 


(130/61) 










Unknown 


sources 








# 


Source 




[Ne 


U] 




[Ne 


111] 


[Fe 


I] 24.04 




src'' 


ext" 


upp" 


src 


ext 


upp 


src ext 


upp 


44 


IRAS 13546-3941 


- 


- 


(27/3) 


- 


- 


(28/4) 


- 


(53/48) 


45 


WL19 


- 


- 


(25/20) 


- 


- 


(19/3) 


- 


(73/30) 


46 


IRAS 16544-1604 














- 


(100/74) 


47 


SSTc2d 1182813.2 


130 


- 


(24/20) 


- 


- 


(26/19) 






48 


SSTc2dJ182849.4 


- 


- 


(100/75) 


- 


- 


(66/57) 






49 


SSTc2dJ182901.8 


- 


- 


(43/6) 


- 


- 


(23/18) 






50 


SSTc2dJ182914.8 


- 


- 


(35/31) 


- 


- 


(29/25) 






51 


EC69 


- 


- 


(4/4) 


- 


- 


(1/1) 


- 


(4/4) 










Disk sources 










# 


Source 




[Ne 


11] 




[Ne 


111] 


[Fe 


I] 24.04 




src-' 


ext" 


upp" 


src 


ext 


upp 


src ext 


upp 


52 


VSSGI 


- 


_ 


(76/49) 


- 


- 


(65/10) 


- 


(130/77) 


53 


IRS 14 


- 


- 


(9/9) 


- 


- 


(6/6) 






54 


OphE-MM3 


- 


- 


(5/5) 


- 


- 


(4/4) 


- 


(63/37) 


55 


GY224 


- 


- 


(28/3) 


- 


- 


(34/5) 


- 


(69/34) 


56 


WL20S 


460 


- 


(44/32) 


- 


- 


(52/8) 


- 


(230/100) 


57 


IRS42 


- 


- 


(150/35) 


- 


- 


(98/17) 


- 


(540/86) 


58 


CRBR 2422.8-3423 


- 


- 


(18/3) 


- 


- 


(16/4) 


- 


(63/47) 


59 


IRS46 


- 


- 


(38/23) 


- 


- 


(36/8) 


- 


(170/78) 


60 


IRS51 


- 


- 


(63/7) 


- 


- 


(55/11) 


- 


(160/66) 


61 


EC82 


- 


- 


(51/9) 


- 


- 


(66/11) 


- 


(110/76) 



src = compact ("or slightly extended, see Table|2t component centered on source ; ext = fully extended ("to Soitzer-IRS) component ; upp = lo" upper limit of src/ext linefluxes 
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Table 3. - continued - 













Embedded sources 












# 




[Fe II] 17.94 






[Fe II] 25.99 


[SI] 






[Si II] 






src 


ext 


upp 


src 


ext upp src 


ext 


upp 


src 


ext 


upp 



I 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 



370 



160 



480 



16 



16 



67 



640 150 



(91/20) 

(50/12) 

(28/25) 

(24/5) 

(19/4) 

(30/20) 

(29/25) 

(100/20) 

(47/45) 

(26/24) 

(25/3) 

(2/2) 

(20/19) 

(110/30) 

(150/38) 

(88/19) 

(16/2) 

(41/38) 

(59/54) 

(36/7) 

(1400/1200) 

(IIO/IIO) 

(2/2) 

(170/150) 

(57/11) 

(88/80) 

(49/41) 

(50/31) 

(85/76) 

(140/130) 

(3/2) 

(4/4) 

(76/60) 

(40/37) 

(10/8) 

(2/2) 

(80/16) 

(6/6) 

(210/190) 

(400/370) 

(160/29) 

(48/36) 

(17/4) 



350 
1300 



360 (83/57) 

920 (300/150) 

(110/75) 



1500 



620 



36 



110 



860 



280 



160 



2700 



260 
150 



(60/40) 

(82/8) 

(93/64) 



(5/5) 

(44/32) 

(130/110) 

(230/150) 

(110/67) 

(48/29) 

(42/22) 

(86/67) 

(89/82) 

(1800/1200) 

(95/7) 

(52/4) 

(210/150) 

(110/93) 

(140/90) 

(57/40) 

(72/59) 

(120/76) 

(130/110) 

(77/69) 

(60/60) 

(160/100) 

(30/20) 

(13/6) 

(80/47) 

(13/13) 

(1400/220) 
(500/340) 



100 



460 
200 



220 



(68/42) 

Unknown sources 



(63/43) 

(150/100) 

(140/83) 



(59/45) 

(61/47) 

(110/58) 



(6/7) 

(45/35) 

(130/94) 

(170/120) 

(110/53) 

(29/18) 

(63/30) 

(78/34) 

(100/92) 

(1600/660) 

(140/78) 

(28/3) 

(260/150) 

(100/83) 

(200/120) 

(74/6) 

(100/67) 

(140/66) 

(150/100) 

(47/47) 

(30/30) 

(240/130) 



670 



160 (44/34) 

260 (580/29) 

(73/64) 

450 (28/27) 

(1700/1300) 
(340/200) 

(77/47) 



2200 



100 



100 



520 



780 



660 



37 
1100 



600 



(99/57) 

(210/140) 

(390/30) 



(190/160) 

(140/91) 

(150/110) 



(39/19) 

(110/82) 

(180/22) 

(340/220) 

(170/130) 

(76/51) 

(110/93) 

(160/91) 

(140/100) 

(2000/1400) 

(170/160) 

(150/13) 

(330/26) 

(240/23) 

(180/140) 

(73/64) 

(93/81) 

(240/170) 

(180/150) 

(470/310) 

(250/250) 

(330/26) 

(67/7) 

(54/63) 

(150/100) 

(56/85) 

(2700/2200) 
(1100/750) 

(150/140) 



[Fe II] 17.94 



[Fe II] 25.99 



[SI] 



upp 



upp 



upp 



[Si II] 



upp 



44 
45 
46 
47 
48 
49 
50 
51 



(38/26) 
(35/34) 

(30/24) 
(44/38) 
(28/25) 
(22/4) 
(11/11) 



(59/39) 

(44/30) 

(120/99) 



21 (3/3) 

Disk sources 



(62/49) 

(54/3) 

(85/70) 



(1/1) 



120 



(67/51) 

(130/12) 

(290/190) 



(12/12) 



[Fe II] 17.94 



[Fe II] 25.99 



[SI] 



upp 



upp 



upp 



[Si II] 



upp 



52 
53 
54 
55 
56 
57 
58 
59 
60 
61 



160 



(69/61) 

(18/18) 

(11/10) 

(40/9) 

(170/39) 

(IIO/IOO) 
(37/33) 
(52/45) 
(77/71) 

(140/130) 



2200 130 



(54/46) 

(30/3) 

(78/6) 

(230/32) 

(140/130) 

(50/41) 

(94/9) 

(75/52) 

(89/73) 



(83/47) 

(35/3) 

(82/73) 

(330/150) 

(120/9) 

(52/4) 

(85/9) 

(140/100) 

(110/85) 



1200 



2000 



(150/140) 

(120/120) 
(110/73) 

(340/230) 
(140/89) 
(77/43) 
(260/29) 
(180/14) 
(160/98) 
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Fred Lahuis et al.: c2d IRS spectra of embedded sources: gas-phase emission lines 



Table 4. Observed linefluxes and Icr uncertainties (10 '^erg cm ^ s ') 



Embedded sources 



# 


Source 


H2 S(0) 


H2S(I) 


H2 S(2) 


H2 S(3) 




src^ ext" 


upp" 


src ext 


upp 


src ext 


upp 


src ext 


upp 


1 


LDN 1448 IRS 1 


_ _ 


(110/71) 


_ _ 


(51/45) 


_ 


(42/37) 


_ 


(140/35) 


2 


LDN1448 NA 


- 


(190/150) 


- 


(31/58) 


77 


(26/3) 


76 


(130/4) 


3 


IRAS 03245+3002 


- 


(240/130) 


14 


(11/2) 


- 


(14/12) 


- 


(13/4) 


4 


L1455 SMMl 






17 


(8/1) 


- 


(14/12) 


14 


(7/3) 


5 


L1455 IRS3 






- 


(18/2) 


15 


(25/3) 


- 


(33/22) 


6 


IRAS 03254+3050 


- 


(80/64) 


- 


(28/23) 


- 


(33/26) 


- 


(52/18) 


7 


IRAS 03271+3013 


- 


(73/58) 


24 


(27/5) 


- 


(32/17) 


25 


(41/3) 


8 


IRAS 03301+3111 


- 


(240/83) 


- 


(100/77) 


- 


(90/72) 


- 


(210/60) 


9 


Bl-a 






- 


(34/5) 


33 


(33/6) 


40 


(38/8) 


10 


Bl-c 






- 


(17/15) 


- 


(20/16) 


- 


(27/9) 


11 


SSTc2d J033327.3 






73 


(13/13) 


17 


(16/2) 


68 16 


(9/2) 


12 


HH211-mm 


150 


(9/8) 


67 


(2/2) 


46 


(4/4) 






13 


IRAS 03439+3233 


- 


(64/39) 


- 


(14/2) 


- 


(14/13) 


68 


(71/13) 


14 


IRAS 03445+3242 


- 


(120/89) 


- 


(67/60) 


- 


(80/45) 


- 


(250/67) 


15 


IRAS 08242-5050 


- 


(250/210) 


- 


(99/80) 


- 


(70/56) 


- 


(810/350) 


16 


IRAS 08261-5100 


- 420 


(130/69) 


- 


(96/61) 


- 


(96/17) 


- 


(200/98) 


17 


Ced 1 10 IRS4 


- 


(54/40) 


23 


(12/2) 


- 


(7/6) 


- 


(37/14) 


18 


Ced 1 10 IRS6 


- 


(76/43) 


- 


(34/30) 


- 


(12/11) 


- 


(210/110) 


19 


IRAS 12553-7651 


- 


(100/100) 


- 


(57/53) 


- 


(40/33) 


- 


(100/17) 


20 


IRAS 15398-3359 


- 


(76/68) 


- 


(23/4) 


- 


(31/29) 


- 


(68/30) 


21 


GSS30-IRS1 


- 


(1700/1100) 


- 


(990/960) 


- 


(540/400) 


- 


(1800/120) 


22 


GY23 


690 


(130/130) 


88 


(74/13) 


- 


(100/85) 


- 


(630/490) 


23 


VLA 1623-243 


66 


(81/7) 


110 


(2/2) 


58 


(3/3) 


- 


(7/3) 


24 


WL12 


- 


(290/190) 


- 


(140/25) 


- 


(130/86) 


- 


(320/190) 


25 


IRS37 


48 


(73/5) 


27 


(33/5) 


- 


(28/23) 


- 


(59/25) 


26 


WL6 


- 530 


(270/91) 


- 


(89/18) 


- 


(110/70) 


- 


(390/65) 


27 


Elias32 


- 


(100/6) 


27 


(36/6) 


- 


(37/34) 


- 


(270/180) 


28 


VSSG17 


62 


(110/8) 


- 


(58/12) 


- 


(88/49) 


- 


(700/430) 


29 


IRS63 


- 


(160/110) 


- 


(65/59) 


- 


(69/39) 


- 


(310/140) 


30 


L1689-IRS5 


- 


(130/120) 


- 


(180/140) 


- 


(100/72) 


- 


(640/240) 


31 


IRAS 16293-2422B 


- 


(99/99) 


66 


(1/1) 


33 


(5/5) 


- 


(4/4) 


32 


IRAS 16293-2422 


110 


(66/20) 


63 


(1/1) 


34 


(1/1) 


- 


(3/3) 


33 


RN0 91 


- 


(230/180) 


- 


(96/87) 


330 


(81/64) 


- 


(530/210) 


34 


SSTc2dJ182916.2 






- 


(39/6) 


- 


(48/36) 


- 


(110/16) 


35 


Seip-S68N 


- 


(49/3) 


44 


(8/9) 


- 


(4/3) 






36 


Seip-SMM4 


^ 120 


(23/10) 


81 


(2/5) 


190 


(5/13) 






37 


EC88 


- 


(110/97) 


- 


(39/12) 


- 


(27/3) 


- 


(4600/4600) 


38 


Serp-SMM3 


97 


(37/10) 


160 


(3/3) 


150 


(3/3) 






39 


R CiA IRS5 






- 


(160/140) 


- 


(150/130) 


- 


(670/260) 


40 


CrA IRS7 A 


- 


(2900/200) 


- 


(600/420) 


- 


(130/120) 


- 


(280/24) 


41 


CrA IRS7 B 


- 


(730/530) 


- 


(120/22) 


66 


(65/13) 


- 


(190/89) 


42 


CrA IRAS32 






- 


(15/3) 


- 


(19/16) 


34 7 


(7/1) 


43 


IRAS 23238+7401 




(80/30) 


Unknown 


(13/8) 

sources 




(10/7) 




(32/9) 


# 


Source 


H2 S(0) 


H2S(I) 


H2 S(2) 


H2 S(3) 




src'^ ext^ 


upp" 


src ext 


upp 


src ext 


upp 


src ext 


upp 


44 


IRAS 13546-3941 


- 


(78/41) 


- 


(31/20) 


_ 


(27/16) 


_ 


(200/140) 


45 


WL19 


21 


(51/3) 


- 


(22/4) 


- 


(20/14) 


- 


(54/18) 


46 


IRAS 16544-1604 


- 


(82/67) 














47 


SSTc2d 1182813.2 






- 


(27/5) 


- 


(26/21) 


- 


(170/12) 


48 


SSTc2dJ182849.4 






- 


(55/46) 


- 


(67/48) 


- 


(630/120) 


49 


SSTc2dJ182901.8 






- 


(35/7) 


- 


(100/99) 


- 330 


(190/69) 


50 


SSTc2dJ182914.8 






- 


(12/11) 


- 


(86/85) 


- 


(130/30) 


51 


EC69 


100 


(3/3) 


60 


(2/2) 


34 


(3/3) 


- 


(20/20) 










Disk sources 










# 


Source 


H2 S(0) 


H2S(1) 


H2 S(2) 


H2 S(3) 




src" ext^ 


upp" 


src ext 


upp 


src ext 


upp 


src ext 


upp 


52 


VSSGI 


750 


(120/110) 


530 


(85/130) 


390 


(88/66) 


_ 


(360/72) 


53 


IRS 14 






790 


(11/11) 


- 450 


(18/18) 


- 


(39/37) 


54 


OphE-MM3 


- 


(41/36) 


3 


(2/0) 


- 


(5/5) 


- 


(lO/I) 


55 


GY224 


30 


(87/5) 


- 


(44/6) 


- 


(29/15) 


- 


(170/25) 


56 


WL20S 


- 


(180/120) 


- 


(76/14) 


- 


(44/37) 






57 


IRS42 


- 


(97/90) 


- 


(68/19) 


- 


(110/63) 


- 


(910/600) 


58 


CRBR 2422.8-3423 


- 


(95/7) 


- 


(27/5) 


14 


(13/3) 


- 


(34/13) 


59 


IRS46 


- 


(130/10) 


- 


(38/7) 


- 


(41/30) 


- 


(250/85) 


60 


IRS51 


- 


(170/130) 


- 


(77/14) 


- 


(35/21) 


- 


(950/580) 


61 


EC82 


- 


(180/140) 


- 


(120/21) 


- 


(52/13) 


- 


(1200/230) 


"see 


footnote a of Table|3J 
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Table 4. - continued - 



Embedded sources 



# 




H2 S(4) 


H2 S(5) 


H2 S(6) 




H2 S(7) 




src 


ext 


upp 


src ext 


upp src 


ext 


upp 


src 


ext 


upp 


1 


- 


- 


(330/40) 


- 


(160/0) 


- 


(160/0) 


- 


- 


(190/0) 


2 


1100 


210 


(250/23) 


440 


(270/57) 


- 


(220/0) 


- 


350 


(670/80) 


3 


230 


- 


(55/11) 


- 


(51/0) 


- 


(100/0) 


- 


- 


(160/0) 


4 


360 


- 


(84/2) 


210 


(34/0) 


- 


(53/0) 


- 


- 


(150/0) 


5 


- 


- 


(100/12) 
















6 


- 


- 


(100/43) 


- 


(320/0) 


- 


(190/0) 


- 


- 


(120/0) 


7 


320 


- 


(48/14) 


550 


(77/0) 230 


- 


(54/0) 


320 


- 


(46/0) 


8 


- 


- 


(320/110) 


- 


(230/0) 


- 


(170/0) 


- 


- 


(190/0) 


9 


910 


- 


(200/15) 


- 


(260/0) 


- 


(140/0) 


3400 


- 


(770/0) 


10 


6000 


250 


(1200/62) 


- 


(480/0) 


- 


(430/0) 


- 


- 


(500/0) 


11 
12 
13 


150 


- 


(10/4) 


- 


(110/0) 


- 


(38/0) 


470 


- 


(45/0) 


_ 


_ 


(90/35) 


_ _ 


(94/0) 


_ 


(100/0) 


_ 


_ 


(90/0) 


14 


- 


- 


(1600/61) 


- 


(510/0) 


- 


(570/0) 


- 


- 


(570/0) 


15 


- 


- 


(950/64) 


- 


(1300/0) 


- 


(1000/0) 


- 


- 


(1100/0) 


16 


- 


- 


(210/68) 


- 


(200/0) 


- 


(210/0) 


- 


- 


(260/0) 


17 


250 


- 


(58/16) 
















18 


- 


130 


(390/25) 
















19 


- 


- 


(1700/100) 


- 


(320/0) 


- 


(410/0) 


- 


- 


(390/0) 


20 


400 


- 


(86/27) 


- 


(160/0) 


- 


(200/0) 


- 


- 


(190/0) 


21 


27000 


- 


(5800/350) 
















22 


- 


- 


(3300/92) 
















23 


- 


- 


(7/4) 
















24 


9800 


- 


(2300/75) 
















25 


- 


- 


(300/35) 
















26 


9000 


- 


(1600/140) 
















27 


- 


- 


(610/23) 
















28 


- 


- 


(2700/35) 
















29 


- 


- 


(2100/54) 
















30 


- 


- 


(1700/290) 
















31 


- 


43 


(5/6) 


67 


(15/15) 


- 


(14/14) 


- 


- 


(15/13) 


32 


- 


18 


(4/4) 


- 


(16/15) 


- 


(15/16) 


- 


- 


(27/27) 


33 


- 


- 


(1300/97) 


- 


(1100/0) 


- 


(1600/0) 


- 


- 


(1600/0) 


34 
35 


- 


- 


(180/130) 


- - 


(150/89) 


- 


(180/120) 


- 


- 


(320/240) 


36 
37 
38 
39 


- 


- 


(8200/8200) 


- 


(250/0) 


- 


(390/0) 


- 


- 


(410/0) 


_ 


_ 


(1300/100) 


_ _ 


(1000/0) 


_ 


(840/0) 


_ 


_ 


(1100/0) 


40 


4900 


- 


(1200/280) 


- 


(3000/0) 


- 


(720/0) 


- 


- 


(690/0) 


41 


4200 


- 


(1000/130) 


- 


(830/0) 


- 


(720/0) 


- 


- 


(790/0) 


42 


90 


- 


(11/2) 


- 


(92/0) 


- 


(60/0) 


- 


- 


(71/0) 


43 


- 


- 


(250/19) 


- 


(130/0) 


- 


(110/0) 


- 


- 


(120/72) 










Unknown sources 












# 




H2 S(4) 


H2 S(5) 


H2 S(6) 




H2 S(7) 




src 


ext 


upp 


src ext 


upp src 


ext 


upp 


src 


ext 


upp 


44 


- 


- 


(290/110) 


- 


(150/0) 


- 


(160/0) 


- 


- 


(170/0) 


45 


- 


- 


(190/26) 
















46 






















47 


- 


- 


(370/70) 


- 


(350/0) 


- 


(410/0) 


- 


- 


(250/0) 


48 


- 


- 


(1100/870) 


- 


(1500/0) 


- 


(2000/0) 


- 


- 


(1500/0) 


49 


- 


- 


(260/110) 


- 


(160/94) 


- 


(260/200) 


- 


- 


(450/380) 


50 


- 


- 


(170/24) 


- 


(310/88) 


- 


(290/100) 


- 


- 


(530/400) 


51 






(35/35) 




(66/65) 

Disk sources 




(64/65) 






(89/86) 


# 




H2 S(4) 


H2 S(5) 


H2 S(6) 




H2 S(7) 




src 


ext 


upp 


src ext 


upp src 


ext 


upp 


src 


ext 


upp 


52 


6600 


- 


(880/120) 
















53 


33 


- 


(8/86) 


- 


(36/0) 


- 


(42/0) 


- 


- 


(37/0) 


54 


- 


- 


(7/6) 


- 


(16/0) 


- 


(22/0) 


- 


- 


(18/0) 


55 


2700 


110 


(190/24) 
















56 






















57 


- 


- 


(2600/62) 
















58 


550 


- 


(120/29) 


- 


(160/0) 


- 


(160/0) 


- 


- 


(210/0) 


59 


- 


- 


(860/44) 
















60 


- 


- 


(1500/36) 
















61 


- 


- 


(910/310) 


- 


(360/0) 


- 


(770/0) 


- 


- 


(520/0) 
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Table 5. Diagnostic parameters of the gas contents 



Embedded sources 



# 


Source 






Warm Hj' 










Hot H2^ 












Compact'' 






Extended 






Compact'' 






Extended 






[K] [1022(,^-2j 


[Mj] 


[K] 


[10™cm-2] 


[Mj] 


[K] 


[10™cm-2] 


[10-^ Mj] 


[K] 


[10'*cm-2] 


[Ur^ Mj] 


1 


LDN 1448 IRS 1 


<600 


< 18.4 




<14 


<7.9 




< 10 


<3.2 




< 2 


< 1.5 


2 


LDN1448NA 


<600 


< 19.2 




<29 


< 17.1 




6 


1.8 


623 


6 


3.8 


3 


IRAS 03245+3002 


<400 


< 12.6 




3 


1.8 




2 


0.9 




0.5 


0.3 


4 


L1455SMM1 


<80 


<2.5 




1 


1.1 




1 


0.4 




0.5 


0.3 


5 


L1455IRS3 


< 164 


<5.1 




<0.8 


<0.4 




<3 


< 1.0 




< 1 


<0.6 


6 


IRAS 03254+3050 


<380 


< 11.7 




<9 


<5.3 




<7 


< 2.2 




< 1 


< 1.1 


7 


IRAS 03271+3013 


<364 


< 11.3 




5 


2.9 


857 


4 


1.2 




0.6 


0.4 


8 


IRAS 03301+3111 


< 1280 


<39.1 




<20 


< 12.0 




< 13 


<4.3 




<6 


<3.6 


9 


Bl-a 


<308 


<9.6 




<2 


< 1.4 




4 


1.5 


615 


3 


1.8 


10 


Bl-c 


< 124 


<3.8 




<6 


<3.8 




6 


2.1 




< 1 


< 1.1 


11 


SSTc2dJ033327.3 


< 104 


<3.2 




10 


5.6 


1514 


1 


0.4 


539 


2 


1.1 


12 


HH211-mm 


<37 


< 1.2 


104 


5 


2.8 




< 1 


<0.4 




3 


2.1 


13 


IRAS 03439+3233 


<244 


<7.5 




<2 


< 1.6 




<3 


< 1.0 




< 1 


<0.6 


14 


IRAS 03445+3242 


<680 


<21.2 




<20 


< 11.1 




<27 


<8.7 




<4 


<2.5 


15 


IRAS 08242-5050 


< 4000 


< 126.5 




<34 


<64.5 




< 152 


<47.1 




< 11 


<21.8 


16 


IRAS 08261-5100 


<2280 


< 70.5 




20 


28.5 




<32 


<9.9 




4 


6.5 


17 


CedllOIRS4 


< 104 


<3.3 




2 


0.8 




1 


0.3 




1 


0.3 


18 


CedllOIRS6 


<212 


<6.5 




<8 


<2.4 




<4 


< 1.3 




<5 


< 1.5 


19 


IRAS 12553-7651 


<300 


<9.3 




< 19 


< 5.5 




<8 


<2.7 




< 2 


<0.8 


20 


IRAS 15398-3359 


<56 


< 1.7 




<5 


<0.5 




1 


0.3 




< 1 


<0.1 


21 


GSS30-IRS1 


<2640 


<81.4 




<309 


<43.4 




48 


14.3 




<25 


<3.5 


22 


GY23 


240 


7.4 




10 


1.4 




< 11 


<3.6 




13 


1.8 


23 


VLA 1623-243 


<29 


<0.9 


142 


<0.9 


<0.1 




<0.1 


<0.0 




7 


1.1 


24 


WL12 


<396 


< 12.3 




< 19 


<2.7 




10 


3.4 




<8 


<1.2 


25 


IRS37 


< 100 


<3.1 




1 


0.3 




< 1 


<0.5 




2 


0.3 


26 


WL6 


<300 


<9.3 




16 


2.3 




10 


3.3 




6 


0.9 


27 


Elias32 


< 120 


<3.7 




0.9 


0.1 




<3 


< 1.1 




5 


0.8 


28 


VSSG17 


< 144 


<4.5 




2 


0.3 




< 10 


<3.3 




9 


1.3 


29 


IRS63 


< 196 


<6.1 




<22 


<3.1 




<7 


<2.3 




<6 


<0.9 


30 


L1689-IRS5 


<288 


<8.9 




<29 


<4.3 




10 


3.2 




<20 


<2.9 


31 


IRAS 16293-2422B 


<26 


<0.8 


> 100 


<8 


< 1.1 




<0.1 


< 0.0 


777 


2 


0.3 


32 


IRAS 16293-2422 


<20 


< 0.6 


109 


<2 


<0.4 




<0.1 


< 0.0 


522 


5 


0.7 


33 


RN0 91 


<356 


< 11.0 




<33 


<5.9 




17 


5.3 




<8 


< 1.5 


34 


SSTc2dJ182916.2 


<400 


< 12.8 




< 1 


< 1.2 




<8 


<2.5 




< 2 


< 1.7 


35 


Serp-S68N 


< 176 


<5.4 




0.7 


0.4 




< 1 


<0.4 




3 


1.9 


36 


Serp-SMM4 


<68 


<2.1 




4 


2.6 




< 1 


<0.5 




11 


6.7 


37 


EC88 


<520 


< 16.4 




<9 


<5.7 




<40 


< 12.5 




<25 


< 15.0 


38 


Serp-SMM3 


<96 


<3.0 


177 


<0.7 


<0.2 




< 1 


<0.3 




34 


7.4 


39 


RCrAIRS5 


<340 


< 10.5 




<65 


<9.9 




< 13 


<4.2 




< 13 


< 1.9 


40 


CrAIRS7A 


< 5200 


< 155.7 




<85 


<22.2 




22 


7.0 




<7 


< 1.9 


41 


CrAIRS7B 


< 1080 


<33.9 




<31 


<8.0 




12 


4.0 




<4 


< 1.2 


42 


CrAIRAS32 


<68 


<2.1 




< 1 


<0.4 


621 


1 


0.5 




<0.3 


<0.1 


43 


IRAS 23238+7401 


<260 


<8.1 




<4 

Unknown 


<2.6 

sources 




<4 


< 1.4 




<0.9 


<0.5 


# 


Source 






Warm H2'' 










Hot H2^ 












Compact'' 






Extended 






Compact'' 






Extended 






[K] [lO^^cm-2] 


[Mj] 


[K] 


[10™cm-2] 


[Ml] 


[K] 


[10™cm-2] 


[10"^ Mj] 


[K] 


[10'*cm-2] 


[10-3 Mj] 


44 


IRAS 13546-3941 


<2440 


<76.0 




<6 


<23.1 




<60 


< 18.2 




<5 


< 19.3 


45 


WL19 


<68 


<2.1 




0.9 


0.1 




< 1 


<0.4 




1 


0.3 


46 


IRAS 16544-1604 


<0.0 


< 0.0 




<0.0 


<0.0 




<0.0 


<0.0 




<0.0 


<0.0 


47 


SSTc2dJ182813.2 


< 184 


<5.7 




< 1 


<0.8 




< 14 


<4.6 




<2 


< 1.3 


48 


SSTc2d J 182849.4 


<88 


<2.7 




< 17 


< 10.3 




< 64 


< 19.9 




< 12 


<7.5 


49 


SSTc2dJ182901.8 


<308 


<9.5 




<0.0 


<0.0 




< 10 


<3.2 




<6 


<3.9 


50 


SSTc2dJ182914.8 


<3 


<0.1 




<4 


<2.6 




< 15 


<4.8 




<2 


< 1.7 


51 


EC69 


<20 


<0.6 


110 


2 


1.7 




<2 


<0.7 




4 


2.4 












Disk sources 














# 


Source 






Warm H2'' 










Hot H2^ 












Compact'' 






Extended 






Compact'' 






Extended 






[K] [1022(,^-2-| 


[Mj] 


[K] 


[10™cm-2] 


[Ml] 


[K] 


[10™cm-2] 


[10-^ Mj] 


[K] 


[lO'^cm-^] 


[Ur^ Mj] 


52 


VSSGl 


284 


8.8 




29 


4.1 




9 


2.9 




15 


2.1 


53 


IRS 14 


<24 


<0.7 


>100 


<90 


< 12.7 




0.5 


0.1 




27 


3.8 


54 


OphE-MM3 


<20 


<0.6 




1 


0.2 




<0.2 


<0.1 




0.3 


0.0 


55 


GY224 


< 120 


<3.7 




1 


0.2 




4 


1.5 




2 


0.3 


56 


WL20S 


<240 


<7.4 




< 11 


< 1.6 




<3 


<0.9 




<5 


<0.8 


57 


IRS42 


< 140 


<4.4 




< 10 


< 1.4 




< 12 


<3.9 




< 14 


<2.0 


58 


CRBR 2422.8-3423 


< 104 


<3.2 




< 1 


<0.2 




1 


0.5 




< 1 


<0.1 


59 


IRS46 


< 140 


<4.4 




< 1 


<0.2 




<4 


< 1.4 




<5 


<0.7 


60 


IRS51 


<220 


<6.9 




< 11 


< 1.6 




<8 


<2.5 




< 10 


< 1.4 


61 


EC82 


< 1200 


<37.7 




< 14 


<8.8 




<35 


< 10.9 




<7 


<4.6 



" 100 K is assumed for the warm component and 1000 K for the hot component. The derived column density and mass of the warm component depend strongly on the assumed temperature. 
A temperature of 150 and 200 K reduces the column density and mass by respectively a factor of ~ 30 and 140. For the hot component 1500 K instead of 1000 K may results in a reduction 
of column density and mass up to a factor of 10. 

''For the unresolved emission a source with r = 50 AU is assumed to obtain an estimate of the column density. 



